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Abstract 



We reconstruct the neutrino mass matrix in the flavor basis, using all available 
experimental data on neutrino oscillations. Majorana nature of neutrinos, normal 
Q_i| mass hierarchy (ordering) and validity of the LMA MSW solution of the solar neu- 

O !* trino problem are assumed. We study dependences of the mass matrix elements, 

m a /3, on the CP violating Dirac, S, and Majorana, p, a, phases, for different values 
of the mixing angle #13 and of the absolute mass scale, mi. The contours of constant 
mass in the p — a plane have been constructed for all m a g. These p — a plots allow 
to systematically scan all possible structures of the mass matrix. We identify regions 
of parameters in which the matrix has (i) a structure with the dominant /ir-block, 
(ii) various hierarchical structures, (hi) flavor alignment, (iv) structures with special 
ordering or equalities of elements, (v) the democratic form. In certain cases the 
matrix can be parameterized by powers of a unique expansion (ordering) parameter 
A w 0.2 — 0.3 (X rd ~ 0.6 — 0.7). Perspectives to further restrict the structure of mass 
matrix in future experiments, in particular in the /3/3o-decay searches, are discussed. 
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1 Introduction 



Significant amount of information about neutrino masses and mixing has already been ob- 
tained from experiments on the atmospheric |IJ and solar neutrinos [|, from laboratory 
experiments, in particular the reactor experiments Q and neutrinoless double beta decay 
searches [||], from astrophysics and cosmology. New substantial results are expected soon. 

What are implications of these results for the fundamental theory, for the mechanism 
of neutrino mass generation, origin of large lepton mixing, relation between the quark 
and the lepton masses? The neutrino masses and mixing appear from diagonalization of 
the neutrino mass matrix. In a sense, the mass matrix unifies the information which is 
contained in the masses and mixing angles (which appear as independent physical observ- 
ables). So, the questions we are asking should be considered in terms of properties of the 
mass matrix. 

It is expected that the structure of the mass matrix can be explained by certain (broken) 
symmetry realized in certain basis at some high mass scale 0. We will call this basis the 
symmetry basis. Thus, to approach the fundamental theory, one should find the mass 
matrix in the symmetry basis and at the corresponding symmetry scale. Both abelian 
(e.g., |7|, U) and non abelian (e.g., 0) symmetries, broken (spontaneously) at the various 
symmetry scales, have been widely considered (see also the reviews ]T(], Also, the 

possibilities have been studied to identify the flavor symmetry scale with other known 
scales in the theory like the Grand Unification scale or the string scale. 

The first step to the fundamental theory is the reconstruction of the matrix in the 
flavor basis using all available experimental data. The flavor basis, formed by v e , v^, v T , is 
determined as the basis in which the mass matrix of charge leptons is diagonal. However, 
the symmetry basis may not coincide with the flavor basis, while the structure of the mass 
matrix depends on the basis substantially. Furthermore, using the existing experimental 
information we can reconstruct the mass matrix at the low (electroweak) scale. The scale 
at which possible flavor symmetry is realized (broken) is unknown. So, the bottom-up 
approach would consist of determination of the structure of the mass matrix at low scale, 
selection of the appropriate symmetry basis and selection of the correct symmetry scale. 

This picture should be taken with some caution: a priory it is not clear whether certain 
symmetry is behind properties of the mass matrix. This should be established by studying 
possible regularities of the mass matrix. 

There is a number of attempts to reconstruct neutrino mass matrix in the flavor basis 



using the available experimental results [fL2| , [13], [lj]. Most of the studies have been per- 
formed in the context of three Majorana neutrinos and the data on atmospheric and solar 
neutrinos as well as from the CHOOZ reactor experiment 0, ^, |4| have been used as an 
input. Clearly this information is not enough to reconstruct the mass matrix completely. 
Apart from the oscillation parameters (mass squared differences and mixing angles), the 
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mass matrix depends on non-oscillation parameters: the absolute mass scale, mi, and the 
CP violating Majorana phases. Furthermore, even not all the oscillation parameters are 
known. In particular, there is only an upper bound on the mixing angle #13 and there 
is no information about the value of the Dirac phase 5. Also the type of mass hierarchy 
(ordering) of the states is unknown. 

The studies performed so far were concentrated, mainly, on identification of the dom- 
inant structures of the mass matrix and possible zeros of certain matrix elements. It was 
realized that in the case of spectrum with normal mass hierarchy, m\ <C m 2 <C m 3 , the 
mass matrix has structure with the dominant /xr-block, formed by M^, M^ T , M TT elements, 



and small elements of the e— row (M ee , M efl , M er ) [T^, [15 1. In the case of inverted mass 
hierarchy, the dominant structure can be formed by elements of the e— row: M e)1 and M er 
P, [12[. These structures may be related to an underlying L e — — L T symmetry. 

In the case of degenerate mass spectrum new dominant structures appear depending 



on the CP-parities of the mass eigenstates (see, e.g., || |l^, |16l). In particular, it has 
been found that the diagonal elements, being equal to each other, can form the dominant 
structure for equal CP-parities of all three neutrinos. Another interesting possibility is the 
dominant structure formed by the ee— , fir— and r/i— elements (moreover, \M ee \ « |M MT |), 
which could imply, e.g., 50(3) flavor symmetry or U(l) symmetry, with charge prescription 
(0,1,-1) and an additional permutation symmetry. Recently, the possibility that some 
matrix elements equal exactly zero has been considered W7\ . 

It was shown that experimental data can be explained in models with universal Yukawa 
couplings |Tj|, which lead to "democratic" mass matrices with all mass matrix elements 
having the same modulus but different phases. 

Completely different approach is based on "Anarchy" of the mass matrix [jnj. It has 
been proposed that the elements of the mass matrix appear as random numbers from 
certain interval and there is no special structure of the mass matrix dictated by certain 
symmetry. It was estimated how frequently neutrino oscillation data can be reproduced 
in this way. Random values of the complex phases of the mass matrix elements M Q/ g have 



also been considered [13 



It was realized that the structure of the mass matrix depends strongly on the unknown 
CP violating phases, especially in the case of degenerate spectrum. In general, in the 
system of three Majorana neutrinos there are three CP violating phases: the Dirac phase, 
S, the unique phase in the mixing matrix relevant for oscillations, and two Majorana phases, 
which are relative phases of the three mass eigenstates. 

In most of previous studies, the CP violating phases were neglected and CP-parities 
have been discussed mainly (see, however, Jl7|] and also [PU| , where the role of phases in 
the generation of large solar mixing is considered). 
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In this paper we perform a systematic and comprehensive study of dependence of the 
neutrino mass matrix structure on the CP violating phases. We concentrate on the first 
step in the "bottom - up" approach: reconstruction of the mass matrix in the flavor basis. 
A short discussion of basis dependence (which deserves a separate study) will be presented 
in sect. 6.7. We suggest a way to analyze all possible structures of the mass matrix which 
are allowed by experimental data. 

The paper is organized as follows. In section |2] we describe our approach and summarize 
physical inputs from neutrino oscillation experiments. In sections |3] and [|, we study the 
dependence of the mass matrix elements on CP violating phases. We consider spectra with 
mass hierarchy (section |3.1| ), partial degeneracy (section f4.2| ) and complete degeneracy 
(section ^^). In section |5] we introduce and describe the (p — a) plots. In section |6| we 
consider implications of CP violating phases for the structure of the mass matrix. In 
section |7] we discuss our result and draw conclusions. 



2 Reconstructing v mass matrix 

The reconstruction of the mass matrix in the flavor basis is the first step of the bottom-up 
approach. The next step - selection of the symmetry basis - requires additional assump- 



tions and therefore is more ambiguous. We will shortly discuss this issue in section pT7 . 
However, already in the flavor basis one can 

- search for regularities in the mass matrix, 

- study correlations of different matrix elements, 

- study correlations between the neutrino mass matrix and charged lepton masses, that 
is study of possible flavor alignment. 

The flavor basis is convenient for searches of symmetries associated with the lepton 
tubers L e , L M , L T . Last but not 
different from the symmetry basis. 



numbers L e , L M , L T . Last but not least, it is not excluded that flavor basis is not much 



2.1 Mass matrix in flavor basis. Parameterization 

The neutrino mass matrix in flavor basis, M, can be written as 

M = U*M dia9 U ] , (1) 

where 

M dia9 = Diag{ mi e- 2ip , m 2 , m 3 e- 2ia ) . (2) 
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Here m« are the moduli of neutrino mass eigenvalues and p and a are the two CP violating 
Majorana phases, varying between and ir. The neutrino mixing matrix U is defined by 



VaL = U ai v iL , a = e,fi,r, i = 1,2,3, 

where v a are the flavor neutrino states, and z/j are the mass eigenstates. We use the 
standard parameterization for U: 

(C13C12 S12C13 s 13 e~ lS \ 

-S12C23 - s 2 3S 13 c 12 e tS c 23 c 12 - s 23 s 13 s 12 e lS s 23 c 13 , (3) 

iS iS I 

S23 s 12 — Sl3 c 23 c 12e — S23 c 12 — <Sl3 s 12C23 e C 23 C\ 3 / 

where = cos 6^, = sin 6^ and 5 is the CP violating Dirac phase. The mixing angles 
vary between and tt/2 and 5 varies between and 27r. 

The matrix M is symmetric and, therefore, defined by six elements^. According to 
Eqs.(|l],@), they can be written explicitly as 

M aQ = 



= {KiU; x ) mie- 2lp + (U* a2 U; 2 ) m 2 + (U* a3 U; 3 ) m 3 e- 2 - , a, {3 = e, /*, r . (4) 

The expression for M^, in terms of rrij, %, 5, p, a, is given in the appendix. 

The mass matrix elements, as functions of CP violating phases, depend on the param- 
eterization. Our choice of parameterization has the following motivations. 

In contrast to previous works, e.g. [T^, 35|, we ascribe the Majorana phases to the 
first and the third mass eigenstates @, so that in the limit of strong mass hierarchy 
the dependence on the phase p disappears. Furthermore, in the limit of degeneracy the 
interplay of the two Majorana phases (due to mixing) is weaker if a is attached to the 
third mass eigenstate. 

We use the standard parameterization (Eq.([3])) of the mixing matrix for two reasons: 
(1) it is the most often used parameterization, in particular in studies of the CP violation 
in neutrino oscillations; (2) the Dirac phase is associated to s 13 . So, the influence of 5 on 
structure of the matrix is suppressed, and moreover, with improvements of bound on s\ 3 
the effect of the phase will decrease. 

Notice that in our parameterization m ee element depends on all three phases. In par- 
ticular, the phase 5 enters in combination a — 5. The dependence of m ee on a — 5 is very 
weak being suppressed by s^ 3 . Another parameterization of the mixing matrix has been 
used, e.g., in fl6| , in which m ee does not depend on S. We find that our parameterization 



Notice that, in general, elements are determined by 6 absolute values and 6 phases. Three phases 
can be eliminated by renormalization of the neutrino wave functions, so that only 9 quantities have 
physical meaning. This matches 9 parameters (3 phases, 3 masses and 3 angles) which appear in our 
parameterization. 
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is more convenient when all elements of the mass matrix (and not only m ee ) are ana- 
lyzed. Moreover, transition to our parameterization in m ee is reduced to the simple shift 
a — > a — 5. 

In what follows we will, mainly, analyze the absolute values of mass matrix elements, 

M a(3 = m^e^ , a, f3 = e, r , (5) 

since the absolute values may give more straightforward information on possible underlying 
symmetry. The phases of the mass matrix elements, (/> a p, are also important for theory. 
The phases are known functions of 9 physical parameters: Q/ g = (p a /3( m i, Oij, P, 8). They 
can be found if these 9 parameters are measured. We comment on the phases and include 
separate discussion in section [ih4| . 

Notice that, in flavor basis, m a p are physical quantities, that is, they can be directly 
measured in physical processes. In particular, the rate of the neutrinoless 2/3-decay is 
proportional to m 2 ee . Other entries are in principle measurable in processes with AL = 2, 
like the decay K + — > 7T~fi + fi + or the scattering e~p — > v e IH^X (for a review see til . 
The rates of these processes are proportional to m^, where a and (3 are the flavors of the 
two produced leptons in the final state or leptons in the initial and final states. 

The present bounds on the elements m a p other than m ee are many orders of magnitude 
weaker than indirect limits. For instance the bound on m 2 efl is 16 orders of magnitude 
above the limit obtained from oscillations 2^] . Clearly, the possibility to improve direct 



limits deserves further studies. 

We introduce the dimensionless quantities 

M af3 = , m af3 = , 

where m 3 is the largest mass eigenvalue. 

2.2 Conservation of the sum of masses squared 

According to Eq. ([]]), the mixing matrix distributes the masses from M dia9 to the elements 
of the flavor mass matrix M: 

m 1 ,m 2 ,m 3 -> m aj3 . 
The following sum rule is useful for analysis of the flavor mass matrix: 

S = m\ = m lp- ( 6 ) 

i=l,2,3 a,(3=e,n,T 

That is, the sum of moduli squared of all the elements of the mass matrix is invariant 
under basis transformation (rotation). The equality (|6]) is the straightforward consequence 
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of the unitarity of transformation. Indeed, denoting M, = M™ 9 (\Mi\ = rrii), we can 
write, using Eqs.([l|J^): 



E m ^ = E 

a, /3 a,P 



J^KiUpM 



E E H<U*pM\ 2 + E E [UaiUpiUajUpjMiM* + h.c] . 
i a,P i>j a,P 

The first term is immediately reduced to £V |Mj| 2 = X^ m i' whereas the second term is 
zero due to orthogonality: U^Uaj = for i 7^ j. 

2.3 Experimental input 

In what follows we will find m Q/ g = m a /3(mi,6ij,5, p,a), using all available neutrino data. 
We will restrict our analysis to the case of normal mass hierarchy (ordering): mi <C 
m 2 <C m 3 {mi < m 2 < m 3 ). The inverted hierarchy (ordering) is disfavored by supernova 



I. Normal hierarchy implies that m 



m\ 



Am^ (m . We will also restrict ourself to the LMA MSW 



SN1987A data @ (see, however, 
and m\ — m\ ~ m\ — m\ 
solution of the solar neutrino problem, which gives the best global fit of the solar neutrino 
data P3| . This solution looks especially plausible after SNO data M and it can be tested 
in the already operating KamLAND experiment |26 ]. We accept interpretation of the 



atmospheric neutrino results [|J in terms of — > v T oscillations as the dominant mode. 
The following experimental information is used. 



The best fit point for the LMA solution [25]: 



tan 2 12 = 0.41, Am 2 sol = 6.2 ■ 10" 5 eV 2 



(7) 



At 99% C.L. the following intervals are allowed: 



tan 2 fl 12 = 0.27 - 0.75, 



(2.7 - 25) • 10" 5 eV 2 



(8) 



From atmospheric neutrino analysis we take, at the 99% C.L. [|TJ: 



tan# 



23 



1 +0.54 
-■--0.35' 



Ami = (2.5±f;l) • 10" 3 eV 2 . 



(9) 



We use the CHOOZ bound on 13 at the 90% C.L. §: 

sin£ 13 <0.22 (Ami > 2 ■ lO^eV 2 ) 



(10) 



7 



In our discussion we will take into account the upper limit on the Majorana neutrino mass 
from neutrino-less 2/3 decay [^], p8| : 



< 0.35 eV (90% G.L.) 



11] 



and m ee < 1 eV, if uncertainties in the nuclear matrix elements are taken into account. 



We think, it is premature to include in the analysis the recent result on 2/5 ^-decay |29] 
which has a controversial interpretation (see discussion in Jl6|, [3(]]). We consider the direct 



kinematic bound on the mass of electron neutrino [pTH , 

m e < 2.2 eV (95% C.L. 



(12) 



The unknown CP violating phases 5,p,a, as well as the absolute mass scale, mi, and 
the angle #13 are treated as free parameters. 

Let us emphasize that the experimental input (0) - ( |10"D does not depend on the Majo- 
rana phases, p and a, and on m 1; because only the differences of raf enter the oscillation 
probabilities. The input does not depend also on the Dirac phase, 5. This can be explicitly 
seen from the parameterization (0): at the level of present experimental accuracy, the so- 
lar and atmospheric neutrino results are determined by U e i, U e 2 and U^, U T ^ respectively. 
CHOOZ gives the bound on \U e ^\. All these quantities do not depend on 5. 



2.4 /ir-block and e-row elements 

In view of large 2-3 mixing, it is convenient to split the six independent elements of the 
mass matrix into two groups: 

- elements of the /ir-block: m^, m Mr , m TT , with zero electron lepton number, L e = 0; 

- elements of the e-row: m efl , m eT with L e — 1, and m ee with L e = 2. 

As we will see later, these groups of elements have different dependences on CP violating 
phases. Moreover, such a split can be motivated by phenomenology. 



2.5 Small parameters, mass ratios and limits 

There are several small parameters in the problem: 
1) The ratio of masses squared differences: 



'A = ,/^-=ai6iss, us) 



atin 



where the central value corresponds to the best fit values of mass squared differences (see 
(0)'®)' an d t ne interval is obtained varying Am 2 atm and Am^, in the ranges given in ([|) 
and (§). 
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2) The 1-3 mixing: s 13 < 0.2 (see ([10])); 

3) The deviation of 2-3 mixing from maximal value (see (^|)), which can be described 

by 

f = cos 29 23 ~ (-0.4 + 0.4). (14) 
Future experimental studies will further restrict all these three parameters. 

Let us introduce dimensionless parameters - the ratios of mass eigenvalues: 

k^, rss^. (15) 
m 2 m 3 

In the case of strong mass hierarchy, k ~ and r = r^. Clearly, we may have k ~ 1 and 
r <C 1. If r ~ 1, then k ~ 1. 

Let us consider the mass matrix in various limits. 

1) r = S13 = £ = 0. In this case we arrive at a matrix with zero e-row elements and 
/ir-block elements equal to m 3 /2. Obviously no dependence on CP-phases appears. 

2) r = S13 = 0, £ ^ (see (Q)-(Q)). We get 

111 

m MM = -(l-0, ^rr = -(l + 0> "V = • ( 16 ) 

The element m AtT is almost unchanged with respect to maximal ^23; while and m rT 
vary with £ significantly and in opposite directions. The determinant of the /ir-block is 
zero. Again, there is no dependence on CP violating phases. 

3) S13 = £ = mi = 0, but r 7^ 0. We have 

2sj 2 r V2s 12 ci 2 r -V2s 12 c 12 r 
M = — I ... e~ 2iCT + c? 2 r e" 2iCT - c 2 2 r I . (17) 



2 



e 2lf7 + c 2 2 r 



Now dependence on the Majorana phase appears in the /ir-block, but there is no phase 
dependence of the e-row elements. 

The influence of the CP violating phases on the matrix structure is very weak in the 
limit of strong mass hierarchy and small S13. Indeed, for r — > 1, the effect of phase a 
disappears, dependence of the elements on the Dirac phase is associated with si 3 , so the 
effect of 5 decreases with S13, the dependence on the phase p is associated to the mass m\ 
and it is negligible when r<l. 
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2.6 Analytic expressions and phase diagrams 



Exact analytic expressions for the mass matrix elements in terms of mass eigenvalues, m^, 
mixing angles and phases are given in the Appendix. We present the matrix elements 
as sums of three contributions corresponding to three different mass eigenvalues in ( |A.l - 
|A.6|) and as series in powers of S13 in ( |A.7|) . Representation of m a p as the sums of three 
terms with different phases is given in (|A.9| , [A.ll| , |A.13| ). We will use various approximate 
expressions for m a/ g which can be obtained from Eqs. (|A.l| - |A.13| ). 

For small S13, one can draw simple graphic representation of the mass matrix elements 
in the complex plane (Fig.|I]). Neglecting terms of order S13 in the brackets of Eqs.( |A.2 - 
|A.6| ) or, equivalently, e terms in Eqs.( [A.9| , |A.11| ), we find that each mass m a p turns out 
to be the sum of three terms with phase factors which depend on certain combinations 
of the phases 5, p, a. So, in the complex plane the masses m a/3 can be represented as 
sums of three vectors (corresponding to the three terms). The lengths of these vectors are 
determined by mass eigenvalues (ratios k and r) and mixing angles. The angles between 
vectors are given by combinations of the phases 5, p, a. 

We will call this graphic representation, used for m ee in |31| and mentioned in |3^, |33 



the phase diagram. The phase diagrams allow one easily to find minimal and maximal 
values as well as phases of the matrix elements, and possible correlations between them. 
In Fig.|I]we show phase diagrams for the case of partial degeneracy: k ~ 1 , r < 1. 

The mass matrix elements are periodic functions of the CP violating phases. In the 
next section we will analyze these dependences by quantifying for each phase the amplitude 
of variations, the period and the average value of the element. The latter we define as the 
average between maximal and minimal possible value of the element. We will also consider 
the relative phases of variations of different elements and correlations between them. 



3 CP phases in the case of hierarchical mass spectrum 



In the limit of strong mass hierarchy, when mi « 0, m\ m Am^, and m\ ~ Am^ (m , k ~ 0, 
only one Majorana phase, a, is relevant. If also s 13 ~ 0, we have, for the matrix of the 
moduli: 

(s\ 2 r c 23 si2C 12 r s 2 3Si 2 ci 2 r 

••• l C 12 C 23 r + S 23 e ~ 2ia \ S23C23 \c 2 12 r - e~ 2ia \ I • (18) 
|„2 e 2 r _i_ „2 2ia\ 
l c 12 6 23' ^ c 23 e I 

Notice that the e-row elements are real, <p ea = 0, whereas for the phases of /ir-block 
elements we have <p a! 3 ~ — 2a. The corrections are proportional to r. 
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3.1 Dependence of m a p on CP violating phases 

In Figs, y we show the six mass matrix elements, rh a p, a, (3 = e, functions of 

the phase a, for different values of the mixing angles 8 23 and #13, from the allowed regions 
given in @ and flUf). Main features of the dependences can be well understood taking the 
lowest order terms in r and Si 3 from Eqs . ( |A . 9| . [A . 1 1| . |A"IT3"D . 



According to (fL8D, the dependence of /ir-block elements on a is a result of the interplay 
of the main, 0(1), term and of the 0(r) term. In the lowest order, the /ir-block elements 
do not depend on S13 (see Figs.||,||); m MT has an opposite phase with respect to the two 
other elements. The relative amplitudes of variations equal 



m 




Am " 

(? l2 r x I tan 2 #23, m TT . (19) 



For the best fit values of the parameters the amplitudes are of order 10%. In the case of 
non-maximal 2-3 mixing, the amplitudes can reach ~ 25%. The corrections ~ rs 13 (|A.1U ) 
lead to small phase shift and small change of the amplitude of variations. 

Neglecting terms of the order rsi 3 s 2 2 , we get from ( |A.11| ) expressions for m eil and m eT : 

m eil w \rs 12 c 12 c 23 + s 13 s 23 e i{S ~ 2a) I , (20) 



m er w |rsi2Ci 2 s 2 3 - Si 3 c 23 e i(<5 2a) \ . (21) 

So, the elements m efl and m er depend on phases in the combination (5 — 2a), they change 
with (S — 2a) in opposite phases, their values are determined by the interplay of the order 
r and order S13 terms, which can have comparable sizes. Maximal values of m eil and m eT 
increase with S13. 

The relative amplitude of variations of m efl with (5 — 2a) is maximal when the two 
terms in ( p0"|) have the same modulus: 

s 13 = s° n = -r sin 2# 12 cot 9 23 w 0.07 . (22) 

If S13 = 5^3, 

m efl = for d — 2a = 7i . 

Maximal value of m eil equals m™ x = 2fn etl = sin2^ 12 c 2 3r. For s 13 < s5 3 , (Fig.||a,fJcJ|c), 
the average value of m eA1 is determined by the first term in (|20"D , whereas the amplitude 
of variations is given by Si 3 /s1 3 . For s± 3 > s° 3 , the second term in ( p0|) dominates. It 
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determines the average value of m eA1 , around which variations occur. The relative amplitude 
of variations is given by the factor s? 3 /si3 (Fig.^JbJ^d,^). 

Behavior of the element m er is similar: the two terms in (21) are equal and can cancel 
each other if 

s i3 = s° 13 = sin 26 12 tan 23 (23) 
(si 3 = Si 3 for maximal 2-3 mixing), so 

m eT = for 5 — 2a = . 

Corrections of order rsi 3 s\ 2l neglected in ( p0[ ) and (|2T|), produce a small relative shift of 
phases of m eAt and rh er (see Fig.0). 
For the ee-element we have: 

m ee ~ \c\ 3 s\ 2 r + sl 3 e 2l{ ^\ . (24) 

It depends on the combination of phases 2(5 — a). Due to the factor rs\ 2 in the first term, 
both contributions in ( p4|) can be comparable in spite of the s^-order of the second term. 
Two terms are equal at tan #13 = S\ 2 ^ « 0.21, that is, near the upper limit for S\ 3 . In 
this case the amplitude of variation can be maximal and 

m ee = for 5 — a = ir/2, 3n/2 . 

Such a situation is approximately realized in Fig.^|b,^|dJ^f. For small values of S13 (si 3 
0.2), the dependence of m ee on phases is negligible (Fig.^]a,^]c,|2]e). The relative amplitude 
of variations is determined by the ratio tan 2 0i 3 /(rs\ 2 ) and the average value equals m ee « 
s\ 2 r. 

Let us analyze the dependence of matrix elements on the Dirac phase 5. The elements 
of //r-block depend on 5 very weakly, via order rsi 3 corrections (see ( |A.10| )). E.g., for 



si 3 w 0.14 (Fig.0b,|d,@f), we find Am^/m^ w Am 5 TT /m TT ~ 0.02. The dependence of 
m^ T on 5 is further suppressed by the factor £ = cos 2# 2 3 . 

The elements of e-row have much stronger relative dependence on 5. As we pointed 
out, the elements m eM and m eT depend on phases in the combination (8 — 2a) (this feature 
is weakly violated by corrections ~ rs 13 , which depend on the phase 5 only). So, up 
to corrections of order s± 3 r, one can extract the information on the 5 dependence of the 
elements from the Fig.[| (or Fig|J for large r) immediately. The change of S by amount AS 
is equivalent to horizontally shift the lines which correspond to m e[1 and m eT along with 
cx-axis by A5/2 and the m ee line by AS, with respect to the lines of /xr-block, which are 
almost unchanged. The phase 5 can be selected in such a way that certain features of the 
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m ee line and other e-row lines will occur at the same value of a. For instance, according 
to Figfjb, one can get m ee -C m efl <C rh eT . 

All the elements have the same period of variation with a, although the phases of 
variations are different. There is a phase shift by n within different groups: 

<j>{m m ) = <p(m TT ) = 0(r?v) + vr = -2a ; (25) 

<K«V) = <j){m eT ) + 7r = -2a + 5 . (26) 

There is a relative shift of phase between /xr-block and e-row elements which is determined 
by 5: 

<K«V) - <K«V/J = 5 , 0(m ee ) - 0(m eM ) = 5. (27) 
These relations are weakly broken by corrections of order rs 13 . 

3.2 Dependence of masses on #12, r and k 

Variations of 612 within the allowed LMA region, given in (§), do not produce substantial 
changes of results shown in Fig.[| With increase of #12, the amplitudes of variations of 
/XT-block elements with a (see ([19])) decrease as c 2 2 . For maximal 1-2 mixing we get 
~ 30% decrease in comparison with the best fit value of 8 12 . In contrast, the amplitude 
of variations of these elements with 5 increases as sin 2^12 (see ( |A.10|) ). The dependence 
on 5 remains weak, because the increase of the amplitude can be only 10%. For the e-row 
elements, the critical value s® 3 is proportional to sin2# 12 (see (p2|) ). The ee-element m ee 
can be two times larger for almost maximal solar mixing angle than for the best fit value 
(see (H)). 

Changes of Am^, and Am 2 atm within the allowed regions, (|8]) and (|9p, produce strong 
effect on the structure of mass matrix. In Fig.[| we show the dependence of mass matrix 
elements on a for r = 0.3, corresponding, e.g., to Am^, w 2 ■ 10~ 4 eV 2 and Am^ m ~ 
2 ■ 10 _3 eV 2 . For the yur-block elements, the amplitudes increase linearly with r (see (P])) 
and for r m 0.3 they can be larger than 30%. For the e-row elements the critical value 
s° 3 (see (p2|)) also increases linearly with r; for r 0.3, we get s° 3 ~ 0.13. For S13 < s" 3 , 
the average values of elements increase as m efl ~ m er ~ r, but the amplitude of variations 
with (5 — 2a) does not change (compare FigfJ, panels a,c,e, with corresponding panels in 
Fig.||D. For S13 > 5^3, the average values of m efJi and m eT do not depend on r, while their 
amplitudes can be maximal (Fig|| panels b,d,f). The average value of ee-element increases 
with r: m ee ~ r; the amplitude of variations with 2(5 — a) does not change. 

Till now, we have considered the case mi = 0. A strong normal hierarchy among mass 
eigenvalues, mi <ti m 2 m 3 , holds for mi up to approximately 0.002 eV (k < 0.3). Notice 
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that, for mi 7^ 0, both Majorana phases become relevant (see (|4])). We have checked that 
varying mi between and 0.002eV, the dependence of m a $ on angles and CP phases, 
showed in Figs|^, is qualitatively the same as for mi = 0, except for the dependence of 
m ee . The ee-element can be about two times larger. Indeed, neglecting terms of order s 2 3 , 
we get: 

m ee m rs 12 (l + k cot 2 6i 2 cos 2p) . 

The second term in the brackets is of order one for, e.g., mi = 0.002 eV, m 2 = 0.006 
eV, tan 2 # 12 = 0.35 and p = 0, it. Depending on p, the ratio of m ee and the other e-row 
elements can significantly change. 



3.3 Structure of the mass matrix in the hierarchical case 

1) As follows from Figs.|2||| the sharp structure with the dominant /ir-block and sub- 
dominant e-row appears for small S13, small r and near maximal 2-3 mixing. In this case 

m RR ~ m ej . ~ m RT < m,, u ~ m,, T ~ m T 



m[e — row) . . . . 

m( l- block) ~M~0.1, (28) 

where m(e — row) and m(/xr — block) refer to typical masses of the e— row and /ir-block 
elements. Improvements of the upper bound on S13 and on £, as well as establishing Am^ 
near its present best fit value, will confirm this structure in assumption of mass hierarchy. 
In the limit of sharp structure [/ir-block]- [e-row], the elements of dominant block depend 
very weakly on S and have about 10% variations (determined by r) due to the phase a. 
The elements m e[1 and m eT depend significantly on the combination (5 — 2a), unless very 
strong upper bound on S13 will be established. The ee— element varies with 2(5 — a), with 
amplitude ~ ^i3' Thus, uncertainties in the structure of the mass matrix due to unknown 
CP violating phases can be substantially reduced by further measurements of mixing angles 
and mass squared differences. 

According to Figs. ||,|3|, for a large part of the parameter space (#23, #13, <5, 0"), the 
structure [/xr-block]- [e-row] is less profound or even disappears. Indeed, in the case of 
large £ or/and large r, the split between masses within the pr-block can be larger than the 
gap between m(e — row) and m(fiT — block), depending on a. Separation of the elements 
in two groups loses any sense. For the extreme case of large values of r, the elements m e[1 
and m er can be even larger than m^ or m TT . 

2) Dependence of the gap between /ir-block and e-row elements on s 13 and r can be 
seen comparing left and right panels in Figs.0,[3] and Fig.0 with Fig.^j, respectively. The 
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deviation of #23 from 45°, leading to a spread among the /xr-block elements (see fll6|)), can 
strongly decrease the gap. 

Let us quantify the size of gap. Taking only leading terms in £, r and s 13 , one has for 
the /ir-block elements: 

m{fiT - block) > m(/ir - block) min = ^(1 - |f | - rc\ 2 ) , 

where m(/rr — block) mm is the value of or m TT , for a = tt/2. The upper bound on the 
e-row elements is given by 

m(e — row) < m(e — row) maa; = —=(sis + rc 12 s 12 ) , 

V2 

where m{e — iow) max is the value of m efl or m eT , for 5 — 2a = or it, respectively. Therefore, 
minimal value of the gap equals 

m(/ir - block)™" - m(e - row) ma * = ^ [l - |£| - r{c\ 2 + V2ci 2 si 2 ) - V2s 13 j . (29) 

One can also characterize the split of the elements by the ratio of mean values of the e- 
row and /ir-block elements. Up to terms quadratic in £, r and S13, m(/rr — block) « rn 3 /2, 
while for the e-row we can take 



/ m % + m er m 3 / 2,2 2 2 \l/2 

m(e - row) = \j — ^ « ^(^13 + r c i2 s i 2 ) 7 

Then 



m e — row / , „ „ „ „ . . . 

^-b.ock) K V 2 W3 + ^), (30) 

in accordance with (^8|) . The ratio ( |30|) does not depend on CP phases and on 623- 

3) Apart from special choice of phases, the ee-element is typically of the order of the 
other e-row elements. 

4) The CP violating phases can change significantly the structure of e-row. As follows 
from Figs. |2],|3|, one can get, e.g., 

m ee <C m efl <C m eT (small r, large Si 3 ) ; 

m ee < m e/1 « m eT (small r, large s 13 ) ; , , 

™ ee ~ "V ~ ™ eT (large r, small si 3 ) ; 

"V < ™ ee < m e r (large r, large si 3 ) . 

Any element of the e-row can be the smallest one. All possible orderings of e-row elements 
can be realized by appropriate choice of the phases. 

5) Depending on phases, one can find a configuration with almost uniform splits among 
the six mass matrix elements and structure with the dominant /ir-block disappears. Still 
the average value of the e-row elements is smaller than the average value of the /ir-block 
elements (see fl3T)|)). Thus one can get flavor alignment (correlation of the neutrino masses 
and masses of charge leptons) . 
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4 CP phases in the case of non-hierarchical 
mass spectrum 

With respect to the hierarchical case, the structure of the mass matrix depends on two 
additional parameters: the mass ration k and the phase p. These parameters enter the 
mass matrix elements in the combinations (see ( |A.7|) ) 

X = s\ 2 ke- 2ip + c\ 2 , Y = s 12 c 12 (l - ke~ 2lp ) , Z = cl 2 ke~ 2ip + s\ 2 (32) 

(in the hierarchical case, k pa , X pa c\ 2 , Y pa Si 2 Ci 2 and Z pa s\ 2 ). We will use the 
following parameterization: 

X = xe i4,x , Y = ye i4>Y , Z = ze i(t>z , (33) 

where x = \X\, <frx, y = |^|, 4>y, z = \Z\, <pz are functions of ^12, k and p. 
In the limit of very small S13 using ( |A.7| ) and the notation (P5|), we have: 

(rz c 23 ry s 23 ry \ 

... \4 3 rx + sl 3 e~ 2i(Tx \ s 23 c 23 \ - rx + e~ 2iax \ , (34) 
\s 2 23 rx + c 2 23 e- 2i<T *\ J 

where ax = <J + <ftx/2- One can first analyze matrix (|34D and then consider corrections 
of the order si 3 . Notice that now the elements of the e-row have non-zero phases which 
depend on p. The dependences of the absolute values and phases of these elements are 
correlated. The phases of /xr-block elements depend mainly on a, with corrections which 
are functions of p. 

4.1 Non-degeneracy case 

For mi < a/ Am^ , we have k < 1, and r <C 1. The largest mass is given by 7723 pa 
a/ Am 2 tm . The contributions of m 1 to the /ir-block elements appear as small corrections, 
but they can be of order 1 for the e-row elements. 

Neglecting terms of order S13, we can use for the /ir-block elements the expressions 
from (0). Comparing with the hierarchical case (see Eq.(|l8|)), we find that the effect of 
mi is reduced to renormalization of the mass ratio r and shift of the phase a: 

x 

r ^r x = r— , a -> a x = + -<Px ■ (35) 
Cf 2 z 

That is, dependence of the elements on phases can be found from Figs.|^J^, by appropriate 
change of r and a. 
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Depending on the phase p, the contribution related to ra\ can suppress or enhance the 
amplitude of variations of pr-block elements with a (see (Jl9|)). The extreme modifications 
are determined by 



rx 



r(l + kt&n 2 612) , for p = 
r(l — ktan 2 6 12 ) , for p = n/2 



For /c < 1 and tan 2 6* 12 < 0.5, the relative effect of mi is below 50%. For p = 0, 7r/2, we have 
0x = and no phase shift occurs. In general, the phase <px is in the interval (— 0^ ax -h0^ ax ), 
where sin0^ ax = A; tan 2 9\ 2 . This maximal phase corresponds to rx = r y/l - k 2 tan 4 9 12 . 
For the elements of e-row, the s 13 corrections should be taken into account (see (|A.7|) ): 



,i{5-2a Y )\ 

m eT « \s 23 ry - s 13 c 23 e l 



"V W I C 23^2/ + Si 3 S23e^ "\ , 

J(S-2a Y )\ 



(36) 



Again, the effect of vii\ is reduced to renormalization of r and a shift of phase (compare 
with Eqs.(H,|2!D): 

r — > ry = r — - — , a — > cry = cr H — 0y . (37) 
•S12C12 2 

Minimal and maximal values of ry are given by: 



ry 



r(l — A;) , for p = 
r(l + fc) , for p = ?r/2 



In these extreme cases there is no phase shift. In general, for arbitrary values of p, <f)y is 
in the interval (— §™ ax -=-(/>y aa: ), where sin^y aa; = k and this maximal value corresponds to 
ry = ryl — k 2 . 

Notice that, for m eil and m eT , modifications of r can be larger than for the elements 
of pr-block; moreover, ry and rx are changing with p in opposite phases. The phases 
of variations of pr-block elements are correlated as in the hierarchical case. No phase 
shift among these elements is induced by mj contribution: in (|25| ) one should substitute 
a — > ax- Similar conclusion is valid for e-row elements: in ([26|) a should be substituted 
by cry. 

For the ee-element, similarly to the previous cases, we get (including s± 3 corrections): 

m ee =\cj 3 s 2 12 r z + s 2 13 e 2i ^ z) \, (38) 

where 



z_ 

'12 



r z = r— , a z = a + -<j) Z . (39) 
S19 z 
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Now the difference between r and r z can be substantially larger: 

f r(l + A;cot 2 12 ) , for p = 
Tz ~ \ r|l - A;cot 2 12 | , for p = tt/2 ' 

and changes with p in phase with rx- Notice that, for k < tan 2 0i 2 , the shift <pz is 
restricted to the interval (— <p r z xax -j- (j)^ ax ), where sin^^"™ = fccot 2 0i 2 . For k > tan 2 i2 , 
the shift is unrestricted. The ee— element is zero for 

, — . . tv 3n 

tan 0i3 = yjrz = s 12 y/r z , (d - a z ) = -, — . 

Since r z can be smaller than r, or even zero, the equality m ee = can be realized for 
smaller values of S13 than in the hierarchical case. Now the strongly hierarchical structure 
of the e-row, 

er j 

can be easily achieved. Maximal value of m ee equals approximately fh^ x ~ r(s\ 2 + kc\ 2 ). 

The dependences of the mass matrix elements on phases can be deduced from Fig.0 
and Fig.|3[ Since, now, the "effective" value of r is different for the pr— block elements 
(rx) and the e-row elements (ry, r z ), one should take, e.g., lines which correspond to the 
pr-block from Fig.|2] and lines which correspond to the e-row from Fig.|3] or vice versa. 

Let us analyze the dependence of the elements on the phase p. The relative amplitudes 
of variations of the pr— block elements are suppressed by a factor s\ 2 rk (see (|A.9|) ): 



cot 2 23 , m m 
h 

1, m, 



* s\ 2 rk x < tan 2 23 , m TT . (40) 

m 

>rr, 



The influence of p on the e— row elements is much stronger. If si 3 ~ 0, we have 



m 



e/i 



C23S 12 c 12 r\l — ke 2ip \ 



(for m eT one should substitute c 23 — > s 23 ) and the relative amplitude of variations is given 
by k. 

The amplitude of ee-element can be maximal if k > tan 2 12 . 



4.2 Partial degeneracy 

For a/ Am^j <mi< ^/ Am 2 tm , we get the spectrum with partial degeneracy m\ « m 2 < 
m 3 . The ratios of masses are 



mi , Am 



2 



< 1 , fc « 1 - . (41) 



AmL ~ 2m 2 
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For mi > 2 • 10 -2 eV, the deviation of k from 1 is smaller than 5% and we can neglect it in 
comparison with other corrections (related to possible large deviations from maximal 2-3 
mixing and to Si% > 0.1). Now the scale of masses is determined by m 3 « a/ m\ + Am^ m ~ 
(1 -T- 2)y / Am^ m . The sum (^|) of all the matrix elements squared equals 

S ~m 2 (l + 2r 2 ). 

Let us consider first the dependence of the masses on phase a (see Fig.f| panels a,c,e 
and phase diagrams in Fig.|l|). In the limit of small S13, we get 

rh TT ^\c 2 23 e- 2i ° + rs 2 23 X 1 (p)\ , (42) 
"V ~ S23C23 \e~ 2la - rX x {p)\ , 

where 

X 1 (p) = = 1) = c 2 2 + 4e- 2i ". (43) 

The mass m MM oscillates with a around s| 3 ; the amplitude of variations depends on the 
phase p. Maximal amplitude is for p = 0, which corresponds to X\ = 1. The mass rh TT 
oscillates in phase with m M(U around the average value c 2 3 ; the mass m MT varies in opposite 
phase. The amplitudes of variations of all /ir-block elements decrease with increase of the 
phase p and it is minimal for p = 7r/2. 

The corrections of order S13 change the amplitudes of variations and produce a phase 
shift. The largest influence of these corrections is for p = n/2. From (|A.9|) we find 

~ l s 23 e " 2iCT + r c| 3 (cos 20i2 - 2e w )| , 
»rr ~ |c^ 3 e~ 2ifT + rs 2 3 (cos 20i2 + 2e rr ) | , (44) 
"V ~ S23C23 |e~ 2ifT - r(cos20i2 + 2e MT )| . 

For 5 = (e > 0), the corrections suppress the amplitude of variations of m M/i and enhance 
the amplitude of rh TT and m Mr variations (see Fig.|]e). For 5 = 71/ 2 (e < 0) the situation is 
opposite: variations of m w are enhanced. 

In the approximation (|42"D, all the elements of /zT-block depend on the phase p in the 
same way. So, there is no relative shift and the relative phases are determined as in (p5|). 
The phase shift seen in FigJ4|c is due to the interplay of e corrections and phase p. 

The dependence of elements of the e-row on a (as well as S) appears due to terms of 
the order S13 (see ( [A. 1 1| ) and Fig.f| panels a,c,e). Neglecting corrections ~ rsi 3 , we get (for 
p being not to close to 0): 

m eil w |rc2 3 Fi(p) + s 13 s 23 e^ 5 - 2 ^ I , m er w |r S 23^(p) - ^se^" 2 ^ I , (45) 
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where 



Yi(p) = Y(k = 1) = s 12 c 12 (1 - e" 2 ^) . (46) 

The masses m eM and m er vary with (5 — 2a) in opposite phase (small phase shift may 
appear due to interplay of order rsi3 corrections and phase p). The amplitude of variations 
is proportional to S13. The average values of the elements increase with p, and they reach 
maxima, m™ ax = rsin2#i 2 c 2 3 and m™ ax = r sin 2#i 2 s 2 3, at p = tt/2. 

Configuration with p = or p m is the special one (see Fig.|]a). In this case the main 
terms in ( |A.11| ) vanish and dependence on phases appears due to e corrections, defined in 



\r sin2# 12 c 23 (e 6jU + e^) - s 13 s 23 e l ( 5 - 2CT ) | 



I \r sin 2# 12 s 23 (e eT + e' eT ) - s 13 c 23 e 1 ^ 2 ^ \ 



(47) 



Notice that elements m eA1 and m eT vary in phase; both the average value and the amplitude 
are proportional to S13. 

Changing the phase 5 by AS one shifts lines which correspond to m efl and rh eT , with 
respect to the lines of /xr-block elements by Aa = A5/2. For instance, according to FigJ4|c, 
one can get equalities m w = m TT = rh^ T and fh etl = rh eT simultaneously. 

Variations of the ee-element (|A.13|) with a as well as with S are strongly suppressed by 
the factor s 2 3 , so that 

m ee «r|Zi(p)|, (48) 

where 

Z l {p) = Z{k = l) = sl 2 + c\ 2 e- 2 ^. (49) 

The average value decreases with increase of the phase p. It varies from m™ ax w r for 
p w 0,7r down to m^ m w rcos2^i 2 for p ps 7r/2. Variations of m ee with p are in opposite 
phase with respect to rh efl and rh eT . 

Let us analyze the dependence of masses on the phase p (Fig.[| panels a,c,e). The 
amplitudes of variations of the pr-block elements with p, AmP oc rsf 2 , are smaller (for 
non-maximal solar mixing) than the amplitudes of o variations. The average values of rh^ 
and m TT decrease whereas the average of m MT increases with increase of a from to ix /2. 
Strong split of masses in the pr-block (see Fig.|5|a,|5|e) is due to cancellation of contributions 
related to m 3 (first term in (|42|)) and to vti\ and m 2 (second term). For large S13, the terms 
of order rs\ 3 can enhance variations with p. 

According to (|45D, variations of the e-row elements with p are strong: the amplitude 
can be close to maximal one. For large values of S13, the phase (2<r — 5) changes signifi- 
cantly the average values of the elements m eA1 and m eT and also modifies the amplitudes 
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of variations with p. 



The matrix elements are all correlated. This can be seen in the limit of very small Sxs- 
For partially degenerate spectrum (k = 1), we get: 



x = z = y 1 — sin 2#i2 sin p 
y = sin 26*12 sin p , 



(50) 



so 



that y = y/1 — x 2 . Moreover, 



, sin2p 

tan0 x = — ■ — . 51 

COt z 0i2 + cos 2p 

The mass matrix can be written as 

rx C23T\ / l — X 2 S23 r \/l — X 2 

rh=\ ... \4 3 rx + s 2 23 e~ 2iax \ S23C23I - rx + e~ 2iax | | . (52) 

S' 2 ry rp _l_ .-i2 ~ 1%(J X I 
1 23' x ' L '23 e 



From (|52[) , we find the following relations among the elements: 

™le + ™% + ™l = r2 ! (53) 

^ = tan ^3; (54) 

m 2 r — m = (1 — r 2 x 2 ) cos 26> 2 3 . (55) 

Notice that m TT = either for #23 = 45° or for rx = 1. The latter corresponds to 
completely degenerate spectrum and p = 0. In this case 

m TT = m^u = \Jl- sin 2 26> 2 3 sin 2 a , m MT = sin 2# 2 3 sin a . 

Furthermore, we find, for the sum of the pr-block elements, 

™ 2 W + »rr + 2rh lr = 1 + T ' ^ \ 

and consequently: 

For a given r the mass matrix (|52[) is determined by ax and x. In general, 0"x and 
x can be treated as two independent parameters. Depending on p, x changes from the 
minimal value x mm = cos 26*12, for p = 7r/2, to x max = 1, for p = 0. The phase </>x varies 
in a rather narrow interval, which decreases with #12: 

sin 4>x ~ ( — tan 2 6*12 -r tan 2 6 ) 12 ) . 
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4.3 Degenerate spectrum 

For mi ^> a/ AmJ tm! we have mi « m2 ~ m3, and the ratio of masses is given by 

^ 2m? " 

For mi = 0.5 eV, the deviation of r from 1 is smaller than 1% and we can neglect it in 
comparison with other small parameters, S13 and £. The e-row elements and the /ir-block 
elements are given by (|A.9| , |A.lT] , [A.13|) , with k = r = 1. Notice that, in the approximation 



Amf (m « 0, the structure of the mass matrix for normal and inverted hierarchy is the 
same. 

Transition to the degeneracy case does not produce qualitative changes in dependences 
of the matrix elements on phases in comparison with partial degeneracy case (see Fig.f| 
b,d,f and Fig.| b,d,f). Amplitudes of variations of the /ir-block elements increase and can 
reach maximal size for specific values of phases. This leads to zero (small) values of certain 
matrix elements and therefore to the appearance of a hierarchical structure of the mass 
matrix. For example, in the case of maximal 2-3 mixing and p = 0, we find from ([52]): 



1 

... l\l + e- 2ia \ i|-l + e" 2iCT | 
i|l + e- 2ifT | 



2 1 



Therefore, m MAt = m TT = for a = it/ 2 (Fig.[|b). Moreover, e corrections cancel in ( A. 9 ) 
for k — 1 and p = 0. For non-maximal 2-3 mixing, using ( f42"D we find that m^ = for 



COS 2#23 _ _ C23 cos 2#i2 + S23 

H • 

-12^ 23^23 



sin p = — — — — — , cos2<r = 2 — 2 

C90 sin 2^19 2cioSo': 



In this case, however, rh TT differs from zero. Such a configuration is realized approximately 
in Fig.|d. 

The average values of m efl and m er increase with respect to the partial degeneracy case, 
whereas the amplitudes of variations with a and p do not change. Average value of the 
ee-element increases with r and can reach 1 for p = (Fig.[|b). 

The amplitudes of variations with p (Fig.^|b,d,f) increase and, for p ps 0, ir, hierarchical 
structure of the mass matrix appears (Fig.|b,|5|f). For some values of phases all the elements 
become approximately equal to each other (see, e.g., Fig.|5]d at p = 1.37r). 

4.4 From hierarchy to degeneracy 

In Fig.|^, we show the dependence of m a p on mi for different values of the Majorana 
phases a and p. As follows from the figure, the hierarchical structure with the dominant 
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/XT-block and small e-row elements exists, independently on phases, for m\j a/ Am^ tm < 0.1 
( m i ^ 0.005 eV). This interval of mi corresponds to hierarchical or non-degenerate spectra. 
The structure with dominant pr— block disappears for m\j y/ Am^ im ~ 0.3 -j- 0.5 (mi < 
(0.02 0.03) eV), that is for partially degenerate spectrum. For mi > a/ Aml tm « 0.05 
eV, the spectrum converges to the degenerate one. In this last case, the structure of the 
mass matrix depends substantially on the Majorana phases. Notice that, in general, the 
pairs of elements ih m and 7Tl TT , 8jS well as m eM and m er , have similar dependences on mi. 

For large part of the phase parameter space, all elements of the mass matrix increase 
with mi being of the same order. Some accidental equalities among them may appear. 
Particular structures are realized for specific values of phases, p, a rs 0, 7r/4, tt/2, shown in 
Fig§ 



5 p — a plots and neutrino-less 2(3 decay 
5.1 p — a plots 

In spite of large freedom related to the unknown CP-phases, a, p, 5, scale mi and S13, 
already the present data give important restrictions on structure of the mass matrix. The 
dependences of various matrix elements on phases are correlated. These features can be 
seen in the p — a plots (Figs.0-|13|) which show contours of constant values of m a p in the 
plane of the Majorana phases p and a. 

Let us comment on properties of the p — a plots. 

The periodicity in p and a implies that the opposite sides of the plots must be identified. 
For example, the case of equal CP parities of Ui, v 2 and 1/3 corresponds to any of the four 
corners of the plots. 

In general, any pair of values (p, a) in the range [0, n) x [0, ir) represents a physically 
independent situation (different mass matrices). However, if 5 = or ir, it follows from (|j) 
that 



m 



P) 



m 



a/3 



V - a, 7r - p) . 



This reflection symmetry is present in Figs.[F||], Figs.|TT|-[T3|, but not in Fig.|TU], where 5 = 
tt/2. 

The phase p is associated with the mass mi, therefore, in the case of strong normal 
hierarchy, the dependence of m a p on p disappears and the iso-mass contours become parallel 
to the axis p. In contrast, the contours for nearly parallel to a axis, since m ee 

depends on a via 0(sl 3 ) terms. There is a relative shift of 7r/2, along the axis a, between 
the patterns for m eM and m eT . 
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The elements m M/1 and m TT have have the same p — a pattern in the limit of maximal 
2-3 mixing and zero S13. The difference between them originates from deviation of #23 from 
45° and from the terms (see ( |A.7| )) 

±sm26 23 s 13 re- i5 Y , (57) 

where the plus sign corresponds to m TT and the minus sign to m^. In the case of maximal 
2-3 mixing, only the term (^) contributes to the difference. The pattern for m^ T is 
complementary to that for and m TT , in the sense that regions of large m^ T correspond 
to regions of small m M/J and m TT and vice versa. 

Small values of the pr-block elements appear at the corners of the plots, p ~ 0, n as 
well as o pa 0, tt, and in the region a ~ tt/2. In the latter case, the corresponding value of p 
depends on 2-3 mixing. For maximal mixing, the regions of small elements are at p ~ 0, ir; 
with deviation from maximal mixing, the regions shift to the center of the plot and merge 
at p ~ Ti / 2 for large values of £. 

Let us comment on specific features of Figs. [7| - [T3| . 

In Fig.[7] we show the plots for the non-degenerate spectrum. There is a sharp separation 
of the e-row and dominant pr-block elements. Structuring within these two groups is rather 
weak. 

In Fig.|8| we show the plots for spectrum with partial degeneracy. Dependence of el- 
ements on p becomes stronger with increase of mi. The pr-block elements have more 
profound structure. The elements m efJi and m er are small in the regions near the corners 
of the plots. 

The plots for spectrum with strong degeneracy are shown in Figs. ||-[T3|. Now the e-row 



elements depend strongly on p, whereas the dependence on a is rather weak. With increase 
of mi the p-dependence becomes stronger for the /ir-block elements (see ( |A.9|) ). The 



patterns for m M/i and m TT differ due to order S13 terms (|57|), which also depend on 5. The 
contribution of the term fl5T|) has minus sign for m MA1 and therefore it adds constructively 
with the other p-dependent term (see ( |A.9| )). For m TT , instead, the contribution has an 
opposite sign, therefore p-dependence remains weak. 

In Fig.|TU| we show the plots for 5 = tt/2. The difference between the plots for m MM and 
m TT becomes smaller in comparison with the case 5 = 0: indeed, for 5 = tt/2, the term 
([57]) has pure imaginary coefficient and its contributions to m MM and m TT become similar. 
For 5 = it, the p — a plots for m MA1 and m TT interchange as compared with those in Fig|| 
The pattern for m^ T is almost unchanged. In the first approximation, the effect of 5 = 71/ 2 
on the e-row elements is reduced to a shift of a by 7r/4 for m efJL and m er and by tt/2 for 



In Fig. [TT] we show the plots for small S13. With decrease of S13, the dependence of 
e-row elements on o disappears, patterns for m w and m TT become more similar, their 
complementarity to the pattern for m^ T becomes sharp. 
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In Fig. [12] we show the plots for non-maximal 2-3 mixing (#23 = 35°). The pattern for 
m ee is unchanged and the one for m MT changes weakly. In contrast, the difference between 
the patterns for m efJi and m er increases. In particular, m efJL can be large for p m tt/2 and 
a ~ 0, it. Also difference of the patterns for and m TT increases. Dependence of m TT on 
phases becomes weaker and regions with very small values of m TT disappear. In contrast, 
for m Pil the region of small values appears near the center of the plot: p ~ a ~ tt/2. For 
6*23 > 45° (not shown) the situation is opposite: region of small values atp~(j~7r/2 
appears for m TT . Also m eT becomes, in general, larger than m efJi . 

In Fig.[13]we show the plots for maximal possible 1-2 mixing. The p dependence becomes 
strong for all the elements and especially for m ee . This element can be zero at p m tt/2. 



5.2 Correlations of mass matrix elements. Extreme values 

The p — a plots allow to systematically scan all possible structures of the mass matrices. 
The pattern of the p — a plots themselves depends on the unknown parameters mi, 5, S13, 
as well as on the uncertainties of the known oscillation parameters. As follows from the 
figures, the dependence of the plots on mi is very strong, whereas the dependences on 5 
and S13 are relatively weak (in view of the strong bound on S13). 

The p — a plots allow to see immediately the correlations between the values of different 
matrix elements. Formally, the 6 independent moduli of matrix elements depend on 5 free 
parameters: 

m Q/3 = m a p(mi,p, a, 5, s 13 ) 

So, only one relation should exist among the matrix elements. Actually, the correlations 
are much stronger due to relatively strong upper bound on s 13 and the fact that the effect 
of 5 is suppressed by a factor s 13 . 

In the physically interesting limits the number of free parameters further decreases. 
Thus, in the case of strong mass hierarchy (mi — > 0) the m± and, consequently, p de- 
pendences disappear: m a p = m a p(a, 5, S13). In the limit of strong mass degeneracy 
the structure of the mass matrix does not depend on the absolute mass scale: m a p = 
mif a p(p,a,5,si 3 ), etc. . 

Each point in the p — a diagram (obviously, the same point should be taken in all six 
panels) corresponds to a mass matrix with certain structure. A given set of p — a diagrams 
(which corresponds to fixed values of m\, S13 and 5) shows 6 elements as functions of two 
parameters: p and a. Therefore, imposing conditions on two (or even one element) one 
may reconstruct whole the matrix up to certain discrete ambiguity. E.g., in the degenerate 
case, imposing condition that m TT is the heaviest element, we find that m ee and m MM should 
be equally large whereas three other elements are small. 

The p — a plots allow to find immediately maximal and minimal values of matrix 
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elements. Using Eq. (|) it is easy to see that the maximal value of the individual matrix 
element is given by 

rn^ = Y,™*\UaiUt>i\. (58) 



It does not depend on the Majorana phases. The minimal value is zero or 

2m&x(mi\U ai U3i\) - y^mi\U ai U3i\ (59) 

i ' ' 

i 

if the latter is above zero. The first term in (|59D is (two times) the largest among the 
contributions from the three mass eigenstates. These statements have been made for m ee 



element in |34| and generalized to other elements in |33], f22 |. 

Extreme values of matrix elements depend on nil as well as on values of oscillations 
parameters. In the case of strong mass hierarchy {mi = 0), the present experimental 
results (see section |3.1| ) allow for zero minimal values of the e-row elements (see also p — o 
plots). Maximal values of the elements are 

m max _ cf 3 S? 2 m 2 _|_ S j 3 m 3 , 

m rnax = C13 [Si 2 (c 23 Ci 2 + S 23 Sl 2 Sl3)™2 + ^IS?^] 

and similarly for m'^ ax , with the substitution S23 *-* c 23 . 

In contrast, the elements of /^r-block have non-zero minimal values. In the limit S13 = 0, 
minimal and maximal values of these elements are, 

and m TT — [s 2 3 <-> c 23 J, — [s 23 , c 23 ^ s 23 c 2 3). 

In the case of strong degeneracy all elements, but m ee , may have zero values. Taking 
into account recent data on solar neutrinos 0, we get m™ n m 0.45m!. Maximal values 
of the /ir-block elements and the ee-element equal: 

m™ ax w m™ x w m™ x w m™ ax w m x . (60) 

In the limit of small S13, maximal values of the two other elements are m™ ax m vii\ sin 2Q\ 2 c 2 3 
and m™ aa; ~ mi sin2^ 12 S23. 

Due to correlations among the mass matrix elements imposed by experimental data 
as well as the sum rule condition (|j), only some elements can take their maximal or 
minimal values simultaneously. In particular, according to the p — a plots of Fig. [7], in the 
hierarchical case only two e-row elements can be zero (very small) simultaneously: m ee 
and m eAt or m ee and m eT . In the case of partial degeneracy m eAt and m er can be very small 
simultaneously. In the case of strong degeneracy, we see, from Figs. |9[-13, that there are 



two groups of elements which can be simultaneously very small: 1) m efl , m eT , m^, m TT \ 
2) m efl , m eT , m^ T . Similarly, from the p — o plots one can get groups of elements which 
reach simultaneously their maxima. 
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5.3 (3(3q v - decay and structure of the mass matrix 

The ee-element, m ee , is the only matrix element for which we have immediate experimental 
access. The p — a plots allow one to find immediately the implications of the results from 
/3/?0i/-decay searches for the structure of the mass matrix (in assumption that the exchange 
of the light Majorana neutrinos is the only mechanism of the decay). For m ee , the Majorana 
phase plots (using a different parameterization) have been considered in |35| . 



The iso-mass contours of nearly parallel to the axis a. Weak dependence of m ee 

on a appears due to term of the order s 2 3 . For very small S13 and (partially) degenerate 
spectrum, the iso-mass contours are determined by 



m ee = m x y 1 — sin 26*12 sin p. (61) 

Suppose that experimental searches give the upper bound m ee < m^. Then, according 
to Figs. |9|-|r5|, there are two iso-mass contours in the p — o plots, which correspond to a 
given value m^ (m ee (p,a) = m^) and a given set of the other parameters (m 1; s 13 , 5, 
etc.): pi = pi(cr) (pi < vr/2) and P2 = P2(c) (pi > 7r /2). The upper experimental limit on 
m ee excludes the following regions in the p — o plots (obviously for all the matrix elements): 

pe[0, Pl (a)] : pG [p 2 (a),7r] . (62) 

The position and the shape of the contours p«(cr) (i = 1,2) depend on mi, 9 12 and S13. 
Taking, e.g., mi = 0.5 eV, S13 = 0.1, tan 2 6*12 = 0.36 and the bound (|iT ), we find from 



Fig.|| that the regions covered by the three darkest strips are excluded. They correspond 
approximately to p < 7r/4 and p > 3ir/4. These regions are excluded for all the elements. 
In this particular case, all corners of the plots and sides with p m 0, n, which correspond 
to hierarchical structure of the mass matrix, are excluded. Clearly no constraint on the 
structure appears for weaker bound, m^ > 0.5 eV (which is allowed by the uncertainty in 
the nuclear matrix elements), or, more in general, for m^ > m\. 

For small S13, the mass matrix can be written immediately in terms of m ee , using 
Eq.©: 



C23 \fr 2 - m? ee s 23 \fr 2 - m? ee 



m=\ ... \c 2 23 m ee + s 2 23 e- 2i ° x \ s 23 c 23 \ - m ee + e" 2 ^ | | . (63) 

I e 2 rn J- r% p~ I 
I 23 6e 23 



Here m ee = xr < r. This form shows how strongly the determination of m ee can influence 
the structure of the mass matrix. The Si 3 corrections to (|63"D, can weakly modify the 
structure of the matrix. 

Positive results of fifto^-decay searches will select two strips in the p — a plot. 
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Substantial bounds on the structure of the mass matrix can be obtained when future so- 
lar neutrino experiments and KamLAND [^] experiment will stronger restrict the allowed 
range for # 12 and also when future (3 decay measurements (KATRIN [P^] ) will strengthen 
the bound on the absolute mass scale. 



5.4 p — a plots for the phases of matrix elements 

The phases of matrix elements are, in general, functions of all the unknown physical pa- 
rameters: Q/ g = (f) a p(mi, p,a,8,sis). In the limits of strong mass hierarchy or/and small 
S13, the expressions are simplified and for some elements the phases are zero. Also in 
certain situations the phases of some elements depend only on p or on a. 

The values of phases correlate (or anticorrelate) with the absolute values of the corre- 
sponding elements. Strong change of phase occurs typically in the regions of parameter 
space where the absolute value of the element is small. There are also correlations between 
phases of different elements. 

In Fig. [14] we show the p — a plots for the phases of matrix elements in the case of 
degenerate spectrum and the same choice of parameters as in Fig. |. Notice that the 
pattern of p — o plots for phases repeats partially the pattern for the absolute values. The 
phase 4> ee depends strongly on p and weakly on a. The phases eM and eT change with p 
and (weaker) with a. The patterns are complementary to some extent: at p ~ 7r/2, <p ef i 
has minimum whereas er maximum. The phases of the /ir-block depend both on a and 
(weaker) on p. The patterns for MM and <p TT are rather similar. Notice that maximal (tt) 
values of these phases are achieved at a ~ tt/2 and minimal (zero) values are at a ~ 0, 7r. 



6 CP phases and structure of the mass matrix 

Possible structures of the mass matrix can be classified in the following way: 

• Hierarchical matrices, with certain dominant and sub-dominant elements. 

• Matrices with flavor alignment 

• Matrices with flavor disorder (flavor "anarchy"). 

• Matrices with certain ordering of elements. In this case, the elements m Q( g have the 
same order of magnitude. 

• Democratic matrices, with equal moduli of all the elements: m a p « mo for any choice 
of a, (3. 

We will discuss these possibilities in order. 
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6.1 Hierarchical mass matrices 



The regions of parameters which correspond to a hierarchical structure of the mass matrix 
can be identified as "white" zones in the p — a plots, where one or several elements have 
small values. Notice that the "white" zones are mainly at the corners or in the center of 
the plot, which corresponds to definite CP-parities or small CP-violating phases. So, the 
most of hierarchical structures can be identified by considering certain CP-parities. 

A systematic search of possible hierarchical structures can be performed in the following 
way. In the limit s\ 3 = 0, the elements of e-row equal: 

M er = -tan6 23 M efl = -rs 23 Y = -rs 23 s 12 c 12 (l - ke~ 2ip ) . (64) 

Since tan 6*23 ~ 0.7 — 1.4, these elements can be either both small or both large. 

Let us consider first the case when m e[1 and m er do not belong to the dominant struc- 
ture, i.e., M efl pa M er « 0. According to fl6~4]), this implies either r — > or p « 0, it. In the 
first case we arrive at the structure with dominant /ir-block: 

/ \ 

M = e" 2w M ... 4, s 23 c 23 \+0(s 13 ,r), (65) 

V c 23 / 

which holds for any value of the phases (see Fig.|7]). Weak ordering of elements is possible 
in the /ir-block. In the second case, p = 0, it, which corresponds to the same CP-parities 
of V\ and u 2 , the ratio r can be of order 1 and new structures appear. For p = 0, tt, we get 
X = Z = 1 and (see (0)): 

/ r \ 

M= ... c 2 23 r + s\ 3 e- 2i ° s 23 c 23 (-r + e- 2i °) + 0(s 13 ) . (66) 
V s 2 23 r + c 2 3 e- 2i ° J 

Such a possibility is realized near the left and right borders of the plots in Fig.[| The 
determinant of the /ir-block is given by 

M TT M^ - Ml w re~ 2 ™ . 

So, with increase of r, it deviates strongly from zero. In the first approximation, we get 
mass matrix with 4 independent dominant elements of the same order: m ee ~ m 



Hierarchy of elements in the /ir-block appears for special values of the phase a. If, e.g. 
tan #23 < 1, we can get M MM w provided that 



a = - , r = tan 2 8 23 
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The mass matrix is then reduced to 

tan 2 #23 \ / r 
-tan #23 ] = (... 
-l+tan 2 # 2; 

Let us underline that such a structure is present in the case of partial degeneracy only. 

In the limit of complete degeneracy, r 1, the condition « requires tan #23 = 1 
and therefore the matrix converges to 

/ 1 0(s l3 ) 0(s 13 ) \ 
rh = ... 0{s\ z ) 1 . (68) 

V OislA I 

This type of matrix has been discussed previously, e.g. in ||. If also 5 = %/2, then Z' = 
(see ( A.8 )) and therefore order S13 terms are zero (see ( |A.7|) and Fig.|lO|). 



If tan #23 > 1, one can get M TT « 0. This, again, requires a = n/2 but r = cot 2 #23 
and, in lowest order in S13, the mass matrix has the form 

cot 2 #23 \ / r 

Mw I ... -l + cot 2 # 23 -cot# 23 = ... -1+r -Vr I . (69) 





It has the same limit (|68]) in the case of completely degenerate spectrum. 

Notice that in the limit S13 — * 0, one should take into account the deviations of k from 
I. This leads to appearance of terms of order Am 2 oZ /2m 2 instead of zeros (see ([ill)). 

According to (|66D , the off-diagonal elements of /ir-block are zero for r = e _2 * CT , that is 



for r = 1 and a = 0, ir. In this case M ee = M w = M rr = 1 and 

/ 1 0( 5l3 ) 0(s 13 ) \ 
m = I . . . 1 (9( S 2 3 ) I • (70) 

So, the dominant structure reduces to the unit matrix (as it has been described, e.g., in 
0) with small off-diagonal corrections (Fig.^|a). If also 5 = 0, it, the order S13 terms are 
zero (see 

Let us study possible equalities of elements of the dominant structure. The conditions 
for the equality = m TT follow from Eq.(|55|). All the elements of /xr-block have the 
same absolute value provided that 2-3 mixing is maximal and a = 7r/4, 37r/4. In this case: 

r 

M= I ... IVTT^e^ IvT + Pe^ |+0(s 13 ), (71) 
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where <fi = arctanl/r. 
For r = 1, we get: 



1 0(s 13 ) 0(s 13 ) 

... 1/V2 1/V2 I . (72) 
1/V2 

Order S13 terms are zero if also 5 = a or a + it. 

Notice that mass matrices considered above depend on r and S23. Dependence on 9 12 
appears only via S13 and Am 2 sol /m\ corrections. 

Let us consider the case where M efJi and M eT belong to the dominant structure. Ac- 



cording to (j64[) , this implies r ~ 1 and p not to close to 0, ir (see regions with p ~ tx/2 



in Figs. (||l3|). In this case, also the element m ee can belong to the dominant structure. 
Indeed, minimal value of m ee is achieved at p = tt/2, so that: 

m ee m ee COS *"l2 



m eil tan 8 23 m eT c 23 sin 26 u 



(73) 



So, all the elements of e-row have comparable values, unless 1-2 mixing is near maximal. 
For 612 ~ 7r/4, the hierarchical structure m ee <C m eA1 ,m er is realized (see Fig.|l3D. 



Let us consider the possibility of zeros in the /ir-block. Now the situation differs from 
that of the case p = 0, it (see Eqs. (|50|j52] )). Since x < 1, we have m Mr 7^ and for maximal 



mixing, all elements of the /ir-block differ from zero. Still, for non maximal 2-3 mixing, 
we can get rh^ = or m TT = 0. For instance, if 823 < 45°, fh^ = when ax = vr/2 and 
tan 2 823 = xr. So, in the case of large e-row elements, one or two of the diagonal elements 
can be zero: m ee , for maximal 1-2 mixing and/or (m TT ), for special relation among 
6>i2, cr and 823 ■ 

Summarizing, the mass matrix has a hierarchical structure: 

(a) In the case of hierarchical mass spectrum: the e-row elements can be about 10 times 
smaller than the /xr-block elements. 

(b) In the case of degenerate mass spectrum: the hierarchy determined by a factor ~ 10 
or more appears in the regions near the corners of the p — a plots: 

P^O-tt/8, (t^O-tt/32, 

for the first corner, and similar interval for the three other corners. In these cases the 
matrix equals approximately to the unit matrix with small off-diagonal terms. Another 
possibility is 

p w - tt/6, a w (0.45 - 0.55) vr 



31 



and similar reflected region p — > ir — p. In this case the mass matrix has a dominant 
structure with m ee pa m^ T , while all other elements are small. 

(c) For non- maximal 2-3 mixing: the element m w or m TT can be small for a pa tt/2 and 
for a value of p which depends on the deviation £ of 2-3 mixing from maximal value. With 
increase of £, the region of small mass approaches the center of p — a plots (p ~ vr/2). 

6.2 Flavor alignment and flavor disorder 

Does matrix show any flavor ordering (alignment), that is, the correlation of the neutrino 
mass terms and the charged lepton masses? To some extent, the lepton mixing matrix itself 
is the measure of the flavor alignment, so that small mixing would imply strong alignment. 
The observed large lepton mixing means weak ordering or absence of the flavor ordering. 
The question of flavor ordering can be studied in terms of mass matrix in flavor basis. In 
this connection, let us consider the possibility that masses decrease with transition from 
the r-flavor to the e-flavor, that is, 

m TT > m MT > m w > m eT > m eil > m ee . (74) 

We will call this possibility the normal flavor ordering or alignment. The ordering with 
m er < iTififi is also possible. Notice that, according to (|55"1), m TT > m wt provided that 
6*23 < 45°. In contrast, one gets from (Q) that m er > m efJi , if 623 > 45°. So, in the 
approximation S13 ~ 0, the "flavor ordering" is impossible. However, for near maximal 2-3 
mixing, the differences (m TT — m^) and (m eT — m efl ) are so small that corrections due to 
non-zero S13 become important. These corrections can produce flavor ordering, as can be 
seen, e.g., in Fig.|3|b, for the case of strong mass hierarchy, mi = 0, and in Fig.|]e (shifting 
e-row lines), for the case r ~ 1. 

There are other possibilities of the flavor ordering. Sets of parameters can be found for 
which the matrix has 

1) rr-alignment (see, e.g., Fig. |3]a, a pa 2.6): 

m TT > m w pa m^ T > m eT pa m efl pa m ee ; (75) 

2) e-alignment (see Fig. |2|d, o pa 2.4), when the masses are sensitive to L e : 

m TT pa m w pa m^ T > m er pa m eM > m ee ; (76) 

3) other alignments, such as (Fig. ||b, a pa 0): 

m TT > m m pa m^ T > m eT pa m eM > m ee . (77) 

Although in many cases m ee can be the heaviest element, inverted flavor alignment 
(when the mass increases with change of flavor from r to e) seems to be impossible. 
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As follows from the Figs. 4 - 6, in a number of cases (partially degenerate, degenerate 
spectrum) the matrix can show flavor disorder. That is, the matrix elements can take 
(relative) values between and 1 without correlation with masses of the charge leptons. 



6.3 Mass matrices with specific ordering of elements 



For mi > a/ Am^ tm (fc « 1), in large part of the phases space, all the elements of the mass 
matrix are of the same order (see FigsJ^|l3|). Values of free parameters can be chosen in 
such a way that any element of the matrix can be the smallest one or the largest one. 
Also one can reach equalities between some of the elements. A number of configurations is 
possible, with only a few restrictions determined by relations among the elements, discussed 
at the end of section p~2| . Varying r, x, ax and #23 (see (p^)), one can get equalities among 
various elements of the matrix. In particular, 



1) m ee = m e[1 for 



C23 

x 



2) m ee = m eT for x given by a similar expression with the substitution C23 <-» S23. 

3) All elements of the e-row are equal for maximal 2-3 mixing and x = l/y/3. 

4) One can reach equality of the diagonal elements m ee = or m ee = m TT and also 



m ee = m m = m TT ; see, e.g., Fig.|c. 



5) The equality of elements of the second diagonal, m eT = m m = m re , is possible, but 
in this case other elements are not small: m TT w m MAt , for example. 

6) According to Fig.[|d, the following equalities can be satisfied: 



m ee = m m = m TT w 2m eM = 2m er = 2m^ T 



for a w 0.7. 

7) For a sa 1.2 (Fig.|d) we find 



m ee = m^r ps 2m TT = 2m eA1 = 2m er = 2m m . 



However, it is not possible to get zero values of all the diagonal elements. Indeed, m ee 
vanishes for r = or x = (the latter corresponds to near maximal 1-2 mixing). However, 
for x = 0, m w and m TT are non-zero: they belong to the dominant structure. The only 
possibility would be to consider inverted hierarchy of the mass eigenvalues. 



6.4 Democratic mass matrix 

It is possible to have equal absolute values for all the matrix elements in the flavor basis. 
To obtain such a "democratic matrix" one should satisfy five equalities among independent 
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matrix elements m a p. In general, we have nine parameters (three masses, three mixing 
angles and three CP violating phases) and we should reproduce the solar as well as at- 
mospheric mass squared differences and mixing angles (4 relations) as well as satisfy the 
CHOOZ bound. So, in principle, the problem is non-trivial. Let us present one realization 
of such a possibility. 

The e-row elements should be as large as the /ir-block elements; this requires r ~ 1 and 
p ~ 7r/2. The /ir-block elements are equal to each other only for a ~ 7r/4, 37r/4. Then, if 
S13 is very small, also £ is required to be very small, otherwise m efJi differs inevitably from 
m eT and the same is true for and m TT . 

Taking the limit S13 = 0, k = 1 (see (0)), we find from equality of the e-row elements: 

£ = , sin 2 29 12 sin 2 p = - . (78) 



According to second equality in (175) the solar mixing angle determines the Majorana phase 



p. Equality of the /ir-block elements leads to the condition 

cos 2<rc 2 2 + cos 2(cr — p)s\ 2 = , (79) 

which fixes the value of a. If conditions ([F_pj) and ( [79| ) are satisfied, it turns out that, taking 
r — 1, the elements of e-row and /ir-block are also equal. Moreover, m af3 = fh = l/y/3. 
Notice that last equality can be immediately obtained from the mass squared conservation 
(section |2T2| ). Indeed, in the degeneracy case we have So/m 2 = ^i^f = 3. For the 
democratic matrix the sum of all elements equals ^ Q ~ m 2 ^ = 9rfiQ. According to the mass 
conservation, we have 9m 2 , = S'o/m 2 = 3, or m 2 , = 1/3. 

6.5 Bi-maximal mixing and its variations 

The Fig.[l3| corresponds to bi-maximal mixing (# 12 = 623 = 45°). Notice that, in contrast 
with pure bi-maximal mixing, 6*13 is non-zero here. The limit #13 — ► leads to disappear- 
ance of dependence of the e-row elements on a and to equality of the patterns for m M/J and 
m TT . 

According to FigJTS], large variety of mass matrix structures can lead to bi-maximal 
mixing. In particular, for p = 0,7r and a = 0, it (corners of the plot), we get the nearly 
diagonal matrix (|70|). For p = 0, ir and a = tt/2, the mass matrix has the form (p5J). For 
p = n/2, it follows x = 0. In this case, neglecting 0(si3) terms, for any value of a we get 
the matrix 



50) 
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discussed in the literature ||. 

Apart from that, many other structures allowed, e.g. matrices with nearly equal ele- 
ments, etc., can lead to bi-maximal mixing. 



Notice that recent data on solar neutrinos strongly disfavor maximal 1-2 mixing 25 



Still mass matrix with bi-maximal mixing can be realized in the symmetry basis. In this 
case the observable non-maximal 1-2 mixing is the result of rotation of the charge lepton 
mixing matrix. 



6.6 Parameterization of M 

Let us consider the possibility to parameterize the mass matrix by powers of a unique 
expansion parameter A< 1: 

m Q /3 = c aP \ n <*0, (81) 

where c a p are numbers of order 1. In the flavor symmetry context, the exponents n a p are 
related to the flavor charges of the corresponding mass terms. If n a p = n a + np, where n a , 
71/3 (a, (3 = e,n,r) are numbers associated with corresponding flavor states, factorization 
occurs: 

m Ql3 = c af3 \ n "\ n e . (82) 
In this case the smallness of various mass terms is correlated: n w = 2n AtT — n TT , 2n efl = 

Hee T^mjii 6tC. 

Let us first consider the case of spectrum with mass hierarchy. As one can see in 
Eq . (|i"8|) , for maximal 2-3 mixing and a w tt/4, 37r/4, all elements of the dominant /ir-block 
can be equal to each other. Then, the elements of the e-row should be suppressed by 
powers of A: 

m^cxA^, p = e,n,r. (83) 

As follows from our analysis, we can have all the e-row elements to be equal among them- 
selves, simultaneously with equality of /ir-block elements: 



moc A 1 1 , (84) 




where (see (|30 



A w \/2(s? 3 + r 2 c\ 2 sl 2) 
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For not too small S13 further structuring is possible, when m ee oc A 2 : 

/A 2 A A \ 

m oc A 1 1 . (85) 




Such a situation is realized, e.g., in Fig|p (with certain shift of the e-row lines due to 5). In 
this case A ~ S13. The matrix (|85|) satisfies the factorization condition. Other structuring 
of the e— row elements is also possible, like (A 2 , A 2 , A) or (A 2 , A, A 2 ) with A ~ 0.3 (see Fig.^j) . 

Mild hierarchy of elements of the fir— block is realized for non-maximal 2-3 mixing 
or/and non-trivial CP phases. According to Fig. |^a, [3]b, we may have m TT « > 
~ frier > fn e/1 m m ee which corresponds to parameterization: 




m oc A 2 A 1 , (86) 



with A ~ 0.3. Also m TT can be the smallest element of the fir— block, instead of m MAt . 

In the case of partial or complete degeneracy, new dominant structures appear and 
therefore new types of expansion is possible. According to Fig.^e and (|71|), the mass 
matrix can have the following form: 



with 



(87) 




A 



«13 



rV2 

Two other possibilities are (see Figs. ||f, ||f): 





A 


M 






A 




(i 


A 




, m oc 




A 2 


■ ) 




1 


Ay 




V A 


1 


A 2 / 



(88) 



with A ~ $13/^/2, which should be taken of order 0.1 for the left matrix and 0.2 for the 
right matrix. 

Notice that value of A which appears in the matrices d84|)-(p8[) and, therefore, consistent 
with present data, can not be too small. We find 

A > 0.1- 0.2 . (89) 
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Values ~ 0.3 — 0.4 are also allowed. The value of the parameter ( |89"D can be equal to sin 9 C , 
where 9 C is the Cabibbo angle, used as an expansion parameter for quark mass matrices. 
In the flavor basis the structure of the charge lepton mass matrix is characterized by the 
two ratios: m^/rrir = 0.059 and m e /m^ = 0.0049. These ratios can also be reproduced as 
powers of A: 

m e : : m T pa A 6 : A 2 : 1 , 

with A w 0.24. 

In a large part of the parameter space, the elements of the mass matrix have the same 
order of magnitude, so that the ratio of matrix elements is close to 1. In this case we can 
introduce the ordering parameter X ord ~ 0(1). Typical value of \ or d can be determined, 
e.g., by the possible spread of /xr-block elements, due to deviation of the 2-3 mixing from 
maximal value: 

X ord pa tan 6> 23 « 0.7 . 

Another possible choice for X ord , in the partial degeneracy case, could be r. We find the 
following structures in FigsfJ,§ (omitting the subscript 'ore?'): 

/ A 4 A 3 A 2 \ / A 6 A 4 A 3 \ 

m oc A 3 A 2 A , m oc A 4 A 2 A . (90) 
\ A 2 A 1 / \ A 3 A I J 

These structures require rather large #13 to enhance the values of the e-row elements. 

In the case of partial or complete degeneracy, situation appears where all elements are 
of the same order with small spread, see, e.g., Fig.^f at a pa 0.7. In this connection one 
can consider the mass matrix as small deviation from the democratic one: 

M oc M D + AM , 

where |Mf^| = 1, AM ~ 0(X) and A is a small parameter. Here A can be taken of order 
S13 or £ or 1 — r (deviation from degeneracy). An interesting possibility could be to take for 
A the deviation of p or a from the values 0, 7r/2, which correspond to definite CP parities. 

6.7 Remarks on the Symmetry basis 

As we have outlined in the introduction, to get further theoretical inference, one needs 
to find the matrix in the symmetry basis and at the symmetry scale. In general, the 
symmetry basis differs from the flavor basis and the mass matrix of charged leptons, Mi, 
is non-diagonal there. The neutrino mass matrix in the symmetry basis, M v , is related to 
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that in flavor basis as M v = UfMUi, where U\ is the mixing matrix which diagonalizes 
M t . 

The matrix JJ\ is unknown and some additional assumptions are needed to fix its struc- 
ture. Clearly this introduces a further ambiguity in the analysis. Here we mention two 
possibilities (two assumptions) which allow one to immediately relate the matrices in fla- 
vor basis and symmetry basis. (The extensive discussion of this issue will be given else- 
where pTp. 

1) It may happen that due to strong hierarchy of the masses of the charged leptons, 
the charged lepton mixing is rather small and U[ « I. In this case, the structures of the 
mass matrix M, discussed in this paper, are not modified significantly under transition to 
the symmetry basis. 

2) Being related to the ratio of masses of the \i and r lepton, the 2-3 angle, Q l 2Z ~ 
y/m M /m T , can be the only large angle in U\ (1-2 and 1-3 mixing angles are very small, 
if they are connected with the tiny electron mass). In this case, effect of charged lepton 
mixing on the neutrino mass matrix is reduced to change of the neutrino 2-3 angle in the 
flavor basis: 

Q _ nsyrn pi 
#23 - £% - (% • 

Taking into account this shift of the angle, one can use neutrino mass matrices obtained in 
this paper as mass matrices in the symmetry basis. This shift can justify large deviations 
of the neutrino 2-3 mixing from maximal value. 

However, there are many models in which charged lepton mass matrix is strongly off- 



diagonal in symmetry basis; see, e.g., [pf|Hl2|p8 



Structures of the mass matrix M will not be modified substantially due to running to 
high scales. It was found that renormalization of M a p is smaller than 10 -4 for the 



Standard Model and about few percents for MSSM. 



7 Discussion and conclusions 

The motivation of our study is to understand how far one can go in construction of the 
theory of neutrino mass using the bottom-up approach, that is, starting from experimental 
results. Neutrino mass matrix in flavor basis unifies information contained in masses 
and mixing angles measured in experiment and therefore can give deeper insight into the 
underlying physics. 

We have elaborated a method which allows one to study dependences of the individual 
matrix elements and of the structure of the mass matrix as whole on the unknown yet 
parameters. In particular, we have performed a systematic and comprehensive study of 
dependences of the neutrino mass matrix elements on the CP violating phases. 
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We have introduced the p — a plots which show contours of constant mass in the plane 
of the Majorana phases p and a. We used the p — a plots to analyze the possible structures 
of the mass matrix. Each point in the p — o plot represents a certain neutrino mass matrix, 
so the p — a plots allow one to scan all possible matrix structures. 

The p — o plots allow to study in rather transparent and straightforward way: 

- influence of the phases on magnitudes of individual matrix elements. In particular, 
one can find ranges in which elements can change and their extremal values (minimal and 
maximal) . 

- correlations between values of different matrix elements. Taking a given element in 
some range one can see immediately intervals in which other elements can change. 

- correlations between the structure of the neutrino mass matrix and the charged lepton 
masses. 

- consequences of experimental measurements of oscillation parameters and m ee on the 
structure of the mass matrix. 

Our results can be summarized in the following way. 

1) The structure of the mass matrix changes significantly with mi. 

For strongly hierarchical mass spectrum (mi « 0) and small S13, the mass matrix has 
a structure with the dominant /ir-block and small e-row elements. The ratio of masses of 
these two groups can be as small as 0.1. 

The dominant structure becomes less profound for large Am^,, large S13 and significant 
deviation from maximal 2-3 mixing. For Am 2 sol > 2 ■ 10 -4 eV 2 , a separation of the elements 
in the dominant /ir-block and sub-dominant e-row has no sense and one can consider certain 
non-hierarchical ordering of the elements. In particular, a configuration with nearly equal 
split among masses is possible. 

For partially degenerate spectrum, the gap between the /ir-block elements and e-row 
elements disappears and all elements can be of the same order. Various equalities between 
the elements and orderings can be realized depending on the CP violating phases. 

In the case of degenerate mass spectrum, the mass matrix can have a hierarchical 
structure with some elements (in particular, from the /ir-block) being much smaller than 
other elements. The hierarchical structures appear for specific ranges of phases. 

In the case of complete degeneracy, the structure of the mass matrix is insensitive to the 
ordering of mass eigenvalues. Therefore, our conclusions are valid also for inverted ordering. 

2) The Majorana phases p and a and the Dirac phase 5 have different impact on the 
structure of mass matrix. This impact depends on values of oscillation parameters and 
mi. 
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(a) The Dirac phase S is associated with the small parameter S13. The influence of this 
phase on the /ir-block elements is relatively weak for any type of spectrum (hierarchical 
or degenerate): it is suppressed by factor s 13 . In contrast, the elements of e-row can 
be substantially influenced by 5, especially in the case of hierarchical spectrum. In the 
first approximation 8 enters m efl and m er in the combination (5 — 2a) and m ee - in the 
combination {25 — 2a). So, the effect of 5 is reduced to the appropriate shifts of phase a 
for rh ee , m ei j, and fh eT . In the p — a plot, for fixed pattern of the /ir-block elements, the 
phase 5 produces a shift of the patterns for m eA1 and m eT , along the axis a. 

Improvements of the upper bound on s 13 in future experiments will further suppress 
the influence of the Dirac phase on the structure of the mass matrix. 

(b) The phase p is associated with the mass eigenvalue mi. So, it has very small effect 
on the mass matrix in the case of hierarchical spectrum. The role of p increases with 
nil. The influence of this phase increases with the solar mixing angle. Therefore future 
measurements of 9 12 in KamLAND and solar neutrino experiments will allow one to further 
restrict the effect of p on the structure of the mass matrix. 

For the best fit value of 6 I 12 , dependence of the /ir-block elements on p is not very strong. 
However, existence of hierarchical structure (zeros) in this block is related to specific values 
of p. There is a strong dependence of the e-row elements on p. Typically m efl and m eT 
have minima at p ~ 0,7r and they are maximal at p pa n/2. The ee-element depends 
on p most strongly. There is a chance to measure/restrict p in the /^/^-decay searches, 
provided that the absolute mass scale will be determined (further restricted) in the direct 
kinematic measurements. 

(c) The phase a is associated with the heaviest mass eigenstate and, consequently, the 
cr-dependence is strong for all the elements but m ee . Variations of the ee-element with a 
are suppressed by a factor sf 3 . 

The phase a enters the e-row elements, m efl and m eT with a factor S13. In spite of this, 
in the case of hierarchical spectrum variations of m efl and m er with a can be strong. With 
increase of r, the relative amplitude of variations of these elements with a decreases. In 
contrast, the dependence of /ir-block elements on a becomes stronger with increase of r. It 
can be enhanced, in addition, if the 2-3 mixing is non-maximal. In the case of degenerate 
spectrum, variations of the /ir-block elements with a can be maximal, so that, at certain 
values of phases, a given element can be zero or the largest one. 

There are correlations among the dependences of the matrix elements on phases. In gen- 
eral, patterns of and 7Tl TT £1X6 complementary to the pattern of m Mr . The patterns for 
m efJi and m eT are shifted by Act = ty/2, etc.. 

3) Using the dependences of the matrix elements on the unknown parameters we have 
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studied possible structures of mass matrices. 

The matrix may have hierarchical form with various dominant structures and small or 
zero elements. The dominant structures can be identified considering the limit S13 — > 0. 
The terms of order S13 give small corrections to the dominant elements. In contrast, the 
Si3-order terms can be important or even give main contribution to the sub-dominant 
elements of the mass matrix. The phase 5 does not determine the dominant structure. 

In the case of hierarchical mass spectrum the dominant structure is formed by the \xr- 
block (see Eq. fl65|) ) . The e-row elements can be about 10 times smaller than the /ir-block 
elements. Properties of this block depend on the 2-3 mixing and on the phase a. 

In the case of degenerate mass spectrum the hierarchy determined by a factor ~ 10 or 
more appears mainly in the left and right-hand sides of the p — o plots. 

One arrives at two rather stable structures: 

(i) the matrix which equals approximately the unit matrix with small off-diagonal 
terms; 

(ii) the matrix which has a dominant structure with m ee ~ m^, while all other elements 
are small. 

Apart from these known hierarchical matrices we have found several new structures 
with non-trivial values of CP violating phases. In particular, for non-maximal 2-3 mixing: 
the element fh m or m TT can be small for a rs it/ 2 and for a value of p which depends on £. 
With increase of £, the region of small mass approaches the center of p — a plots (p ~ 7r /2). 

Typically, CP violating phases which differ substantially from from , ir/2 or 7r lead 
to non-hierarchical matrices. 

We have found that the matrix may have certain flavor ordering (alignment), when 
masses increase with change of the flavor from e to r. At the same time we find that the 
data can be reproduced by matrices with flavor disorder, when no correlation between the 
size of the mass terms and the flavor is observed. The democratic mass matrix is also 
possible. 

4) Typical separations among the elements in the hierarchical structures of the neutrino 
mass matrix are characterized by a factor 0.2 - 0.3. We have found that it is possible to 
parameterize the matrix by powers of a single parameter A (whose origin can be in the 
breaking of some flavor symmetry at high energy). The value A ~ 0.2 — 0.3 is consistent 
with the Cabibbo angle and also it can be related to the ratios of charge lepton masses. 

If 2-3 mixing is not maximal, one can introduce an ordering parameter \ or d ~ tan #23 ~ 
0.6 — 0.7. We find that the whole matrix can be parametrized in terms of powers of this 
ordering parameter. 

5) The following results from forthcoming experiments will have crucial impact on the 
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structure of the neutrino mass matrix: 

- improvement of bound on (or determination of) the deviation £ from maximal 2-3 
mixing; 

- precise determination of the solar oscillation parameters, Am^ and 612] 

- improvement of bound on (or determination of) S13; 

- improvement of bound on (or determination of) m ee ; 

- direct kinematic measurements of the neutrino mass. 

Is it possible to determine uniquely the mass matrix, at least in principle? The answer 
depends on future experimental results. Let us take the most optimistic situation: suppose 
that the neutrinoless 2/3 decay is discovered with m ee > 0.1 eV and the direct measurements 
of neutrino mass give m > 0.5 eV with high precision. Let us assume also that mixing 
angles are measured with high accuracy. In this case, the spectrum is strongly degenerate 
and one can use neutrinoless 2j3 decay data to determine the CP violating phases. The 
problem is that m ee depends both on p and on a, and moreover, the dependence on a is 
very weak being suppressed by sf 3 . This means that p can be measured with rather good 
accuracy, whereas no bound on o can be obtained: small variations of p can imitate effect 
of a in the whole possible range. The only exception is if the measured m ee is at maximal 
(or minimal) possible value predicted for a given (measured) absolute scale of the mass. 
That would correspond to certain CP-parity of v\ and 5 — a = or tt/2. Then, measuring 
5 (in neutrino oscillation experiments) one can get o. Clearly even this program looks very 
challenging. Other experimental situations are even more difficult. 

The determination of a looks practically impossible, unless methods of direct measure- 
ment or independent reconstruction of at least one another matrix element (apart from 
m ee ) will be found, 

The p — a plots give an idea of uncertainty in the structure of the mass matrix if a is 
unknown. If 2(3 decay searches give a positive result and direct measurements improve the 
bound on (or measure) m 1; we will be able to select a narrow vertical strip in the p — o 
diagram. This will also restrict other elements, but significant uncertainty will be left due 
to their dependence on the phase a. In particular, as follows from the figures in the case 
of degenerate spectrum, the structure of the /xr-block will be largely unfixed. 

The hope is that even partial reconstruction of the mass matrix may give important 
hint in favor of certain underlying theory. 
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Appendix: Analytic expressions for matrix elements 



For convenience of the reader we present here explicit expressions for the neutrino mass 
matrix elements in flavor basis, M a p = e,/i, r), as functions of the mass eigenvalues 
rrii , of the mixing angles % and of the CP violating phases <5, p, a (see (§)). 

Denoting = cos 6 \j and Sij = sin 6^-, we get the following expressions for the matrix 
elements, as sums of three terms corresponding to the three mass eigenstates: 



c i3 c\ 2 mie 2%p + 

C 13 S 12 m 2 + 

s 2 l3 m 3 e 2i ^- a) ; 



(a.i; 



M, 



Cl3 C12 (-C 23 S12 - S 23 C U S 13 e~ lS ) Td\ e~ 2tp + 

ci 3 si 2 (c 23 C12 - s 23 S12 s i3 e~ J(5 ) m 2 + 
C13 s 2 3 si 3 m 3 e i( - s ~ 2 ^ ; 



(A.2) 



M er = c i3 cia (s 23 S12 - c 23 c i2 s 13 e _J<5 ) m 1 e~ 2lp + 
ci 3 S12 (-s 23 ci 2 - c 23 s i2 s i3 e _4<5 ) m 2 + 
ci 3 c 23 si3 "23 e i{5_2<T) ; 



(A.3) 



M, 



i'i-i 



c 2 3 S12 - s 23 ci2 si 3 e l<5 ) mie" 
S23 S12 S13 e" J<5 ) m 2 + 



2ip _|_ 



(C23 c 12 



c? 3 4s m 3 e 2iCT ; 



Af 



(s 23 S12 - c 23 C12 S13 e * 5 ) (-c 23 si 2 - s 2 3Ci2Si 3 e tS ) m 1 e 2lp + 
(-S23 C12 - c 23 S12 s 13 e~ J(5 ) (c 23 C12 - s 23 S12 si 3 e~ 4<5 ) m 2 + 



c? 3 c 23 S23 ^3 e 2i<T ; 



(A.4) 



(A.5) 



M TT = (s 23 sia - c 23 C12 S13 e- i<5 ) 2 m 1 e~ 2ip 
(-S23 C12 - c 23 S12 si 3 e~' i<5 ) " m 2 + 
c 2 3 c^ 3 m 3 e~ 2ia . 



(A.6) 
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It is convenient to use also a representation of the matrix elements as series in powers 
of S13. Using the equations ( |A.1| )-( [A~6| ) and the definitions (0), we get: 



M ee = rZ + s\ 3 Z> , 

M eil = c 13 (c 23 rY + s 13 s 23 e- iS Z') , 

M eT = c 13 {-s 23 rY + s 13 c 23 e- i5 Z') , 

= c\ 3 Xr + s 2 23 e~ 2i(T - sin29 23 s 13 re- i5 Y - s 2 23 s\ 3 e- 2i& Z' , 
M TT = s 2 23 Xr + c\ 3 e- 2ia + sin 26 23 s 13 re~ iS Y - c 2 23 s 2 13 e~ 2i6 Z' , 

= s 23 c 23 (-rX + e" 2iCT ) - cos 26 23 s 13 r e - i5 Y - s 23 c 23 s 2 3 e- 2lS Z' 



(A.7) 



where 



Z' = e 2i{5 ~ a) - rZ . (A.8) 



We will use also another form for the matrix elements, which is obtained neglecting 
terms of the order s 2 3 . For the /XT-block, we get: 

TTlfifi ~ \ s 23 e rC 23( C 12 ~ ^W*) ^ rC 23(' S 12 ^/J./J-) e P \ i 

4 \c 2 3 e- 2l ° + rs 2 23 (c 2 2 + e TT ) + krs 2 3 (s 2 12 - e TT ) e - 2ip \ , (A.9) 



m 

rh^ » S23C23 |-e" 2iff + r(c 2 2 + e^) + fcr(s? 2 - e^ T )e- 2ip \ 



where 



tan #23, 

e aP {5) = sin29 12 s 13 e~ i5 x <j cot 023, • (A.IO) 

COt 26*23, V 



Notice that the three different terms in the Eqs.( |A.9| ) are contributions of the three masses 



777,3, m 2 and mi. For the elements of the e-row we get, up to an overall factor C13: 



m P 



|rc 12 s 12 C23 [1 - e 6(U - he'^l + e^)] + sia^e^" 2 ^ | , , 
\rc 12 s l2 s 23 [l + e eT -ke~ 2i P{l-e' eT )}-s 13 c 23 e^~ 2 ^\, 1 ' ' 



where 



(e efl , e' e/ J = s 13 e t5 t&n9 23 x (tan6 12 , cot9 12 ) 



(A.12) 



e' eT ) = si 3 e lS cot #23 x (tan0i 2 , cot0i 2 ) • 
The ee-element can be written as 

m ee = \c\ 3 (c\ 2 ke- 2ip + s 2 12 )r + s\ 3 e 2l{5 - a) \ . (A.13) 
It does not depend on the angle 6*23. 
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Figure 1: Phase diagrams. Shown is a graphic representation of the mass matrix elements m a /3 
in the complex plane. Thick lines represent the contributions of the three leading terms in 
expressions C |A79| , |A. 1 HjA7l3|) . The dia gram for m TT is obtained from the m w diagram with the 
substitution S23 <-> C23. The length of the dash-dotted line gives the value of the matrix element. 
The diagrams correspond to the spectrum with partial degeneracy (k ~ 1, r < 1). 
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Figure 2: Dependence of the absolute value of neutrino mass matrix elements (in units 
\J Am^ tm = 0.05 eV) on a, for different values of 623 and #13. We take tan 2 9±2 = 0.36, 6 = tt/3, 
m,2 = 0.14 7B3, mi = 0. 
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Figure 4: Dependence of the absolute value of mass matrix elements (in units \J Am 2 im ) on 

a, for partially degenerate spectrum (panels a,c,e) and completely degenerate spectrum (panels 

b, d,f). We show dependences for different values of the phase p. We take Am^ = 5 • 10 _5 eV 2 , 
Am 2 atm = 2.5 • 10~ 3 eV 2 and tan 2 12 = 0.36, tan6» 23 = 0.93, s 13 = 0.1, 5 = 0. 
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a: cr=0, m 1= 0.03 eV b: a=0, mi=0.5 eV 




Figure 5: Dependence of the absolute value of mass matrix elements (in units \J Am^ m ) on 
p, for partially degenerate spectrum (panels a,c,e) and completely degenerate spectrum (panels 
b,d,f). We show dependences for different values of the phase a. We take Am^, = 5 • 10 _5 eV 2 , 
Am 2 atm = 2.5 • l(T 3 eV 2 and tan 2 n = 0.36, tan<9 23 = 0.93, s 13 = 0.1, 5 = 0. 
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Figure 6: Dependence of the absolute value of neutrino mass matrix elements (in units \J Am 2 tm ) 
on m\. We show dependences for different values of the phases p and a. We take Am^ = 
5 • 10~ 5 eV 2 , Am 2 atm = 2.5 • l(T 3 eV 2 and tan 2 9 12 = 0.36, tan6> 23 = 0.93, s n = 0.1, <5 = 0. 
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Figure 7: The p — a plots for non-degenerate spectrum, with m\ = 0.005 eV. Shown are contours 
of constant mass (iso-mass) m = (0.1, 0.2, . . . , 0.9)m max , where rn max = 0.03 eV is the maximal 
value that the matrix elements can have, so that the white regions correspond to the mass 
interval (0 - 0.003) eV and the darkest ones to (0.027 - 0.030) eV. We take Am 2 sol = 5 • 10" 5 eV 2 , 
Am 2 atm = 2.5 • 10" 3 eV 2 and tan 2 12 = 0.36, tan6> 23 = 1, s 13 = 0.1, 5 = 0. 
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Figure 8: The p — a plots for partially degenerate spectrum, with mi = 0.03 eV. Shown are 
contours of constant mass (iso-mass) m = (0.1, 0.2, . . . , 0.9)m max , where rn max = 0.045 eV, so 
that the white regions correspond to the mass interval (0 — 0.0045) eV and the darkest ones to 
(0.0405 - 0.045) eV. We take Am 2 sol = 5 • 10~ 5 eV 2 , Am 2 atm = 2.5 • 10" 3 eV 2 and tan 2 9 12 = 0.36, 
tan #23 = 1, si3 = 0.1, 5 = 0. 
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Figure 9: The p — a plots for completely degenerate spectrum, with mi = 0.5 eV. Shown are 
contours of constant mass (iso-mass) m = (0.1, 0.2, 0.9)m max , where rn max = 0.5 eV, so 
that the white regions correspond to the mass interval (0 — 0.05) eV and the darkest ones to 
(0.45 - 0.5) eV. We take Am 2 sol = 5 • 10~ 5 eV 2 , Am 2 atm = 2.5 • 10~ 3 eV 2 and tan 2 12 = 0.36, 
tan #23 = 1, S13 = 0.1, 5 = 0. 
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Figure 10: The same as in Fig.||, but for non-zero Dirac phase: 5 = tt/2. 
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Figure 12: The same as in Fig.^, but for non-maximal 2-3 mixing: 623 = 35°. 
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Figure 14: The p — a plots of the phases of mass matrix elements in the case of degenerate 
spectrum. We use the same set of parameters as in Fig.^. Shown are the contours of constant 
phase (iso-phase) cp a ^ = — ir + mr/8 for n = 0, 7 (dashed contours) and for n = 8, ...,15 
(continuous contours). 



62 



